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INTRODUCTION 


+ ee AERODYNAMICS of an airplane includes the 
consideration of all factors directly affeeting the 
performance, stability, and control of the airplane. 
Since stability and control qualities are designed to 
meet certain requirements for satisfactory flight and 
are then fixed in the airplane, these qualities are not 
included in the term ‘‘operational aerodynamics,” 
whieh implies a choice of the manner in which an 
airplane’s characteristics can be used to obtain the 
most efficient operations. Operational aerodynamics 
shall then be taken to mean the study of the optimum 
use of an airplane’s performance abilities. 

The words ‘“‘optimum”’ and “efficient,’’ when applied 
to an airplane’s operation, ordinarily imply obtaining 
the lowest cost per pound-mile. However, other 
characteristics must be qualitatively considered. In 
addition to the actual effect of speed on cost, there is 
an indirect effect of speed in terms of traffic appeal. 
Ride roughness is another semi-intangible quality, 
which is important but which can be treated only by 
establishing acceptable limits. 
operation used may have an influence on air traffic 
control, which has a direct effect on cost and speed 
but which would be difficult to evaluate in specific 
terms. Therefore, it is generally necessary to study 
the operation of an airplane in terms of the quantities 
whose costs we can estimate and then attempt to select 
the most practically efficient operation by balancing 
the unmeasurable advantages and disadvantages of 
various methods of operation against the calculable 
costs: 

For transport planes in general use until recently 
the pptimum use of the airplane was easily determined. 
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Using the DC-3 as an example, the maximum cruising 
altitude is limited to about 14,000 ft. by the lack 
of a supercharged cabin, so that only moderate altitude 
variations are possible. The rate of descent is limited 
by the maximum rate of increase of air pressure ac- 
ceptable for passenger comfort, to about 400 ft. per 
min. This limited rate of descent automatically re- 
stricts the air speed in the descent and keeps the air 
speed well under the placard glide speed, 257 m.p.h., 
indicated. 

Modern airplanes in the DC-6 class have a wider 
variation of operation possibilities. The passenger- 
comfort ceiling is greatly raised, since at 23,000 it. 
the pressure in the supercharged cabin is equal to that 
at an atmospheric altitude of 10,090 ft. Higher 
engine critical altitudes make possible improved per- 
formance at high altitudes. The level flight speeds 
have so increased that the structure is designed for a 
placard glide speed of 360 m.p.h., indicated, thus con- 
siderably increasing the speed structurally possible in 
a descent. The rate of climb and rate of descent per- 
missible by passenger-comfort considerations are in- 
creased because of the pressure cabin, the exact rate 
being dependent on the cruising altitude. Selection 
of optimum methods of climbing and descending is 
sufficiently complex for this airplane to require a study 
of the effects of the operational parameters on block 
speed, operating costs, and ride roughness. 

The next step will be toward gas turbine powered 
aircraft, either jet or propeller driven. The trends 
seen in the DC-6 as compared to the DC-3 are now 
greatly accentuated. The supercharging of the cabin 
will be greater, since high altitudes are necessary for 
gas turbine flight for satisfactory fuel economy. The 
altitude range will extend up toward 30,000 to 40,000 
ft. The greater degree of cabin supercharging will 
increase permissible rates of ascent and descent. If 
the high level flight speeds that are now possible are 
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PLAN VIEWS 
O0C-3, DC-6 AND JET TRANSPORT AIRPLANES 


JET_TRANSPORT 
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to be permitted at low altitudes, the airplane will 
have to be stressed for considerably higher indicated 
speeds. The possible cruising altitudes and rates of 
climb and descent now become so varied that a de- 
tailed study is clearly necessary to determine the most 
efficient method of operation compatible with safety 
and comfort requirements. 


METHOD oF STUDY 


Calculations of the range, cruising speed, block speed, 
and operating costs have been carried out for various 
extreme conditions of cruising power, altitude, rate of 
climb, and rate of descent operation. The effect of 
the operational variables on the cruising characteris- 
tics have been studied for the DC-3, the DC-6, and a 
typical large turbojet airplane. Plan views of the 
three airplanes are shown in Fig. 1. The DC-3 has 
a maximum take-off gross weight of 25,200 Ibs., has 
a wing area of 987 sq.ft., and is powered by two engines 
with a maximum continuous rating of 1,050 b.hp. 
each. The DC-6 has a maximum take-off gross weight 


of 93,200 Ibs., although for the particular operation | 


studied here the take-off weight is about 84,000 Ibs. 
It has a wing area of 1,462 sq.ft. and is powered by 
four Pratt & Whitney R-2800-CA 15 engines with a 
maximum continuous rating of 1800 b.hp. per engine 
in low blower and 1,600 b.hp. per engine in high blower. 
The jet-powered transport has a design take-off weight 
of approximately 76,000 Ibs., has a wing area of 1,100 
sq.ft., and is powered by four turbojet engines with a 
rated maximum contintious static thrust at sea level 
of 4,960 Ibs. per engine. The relative effects of the 
operational variables on the three airplanes have been 
compared and are presented in graphical and tabular 
form. The methods used are summarized below. 


Rate of Climb 


‘ Rate of climb has been computed using maximum 
continuous power and a maximum cruising power for 
the DC-3 and the DC-6, and maximum continuous 
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and a maximum cruising thrust for the jet transport, 
The maximum continuous power and thrust are 
called rated power and thrust, respectively, in the dis- 
cussion in this paper to clearly differentiate between 
these and maximum cruising power and thrust. Since 
a preliminary study of the effect of climbing speed 
had indicated that optimum block speed and economy 
were obtained at about the speed for best climb, this 
climb speed was used for comparison of the three types 
of airplanes. No limit on rate of climb from a pas- 
senger comfort standpoint was established, since studies 
have shown that cabin pressure rates of climb up to 
1,000 ft. per min. are not objectionable and none of the 
airplanes studied exceeded this value. 


Cruising Operation 


Cruising has been assumed to be at the chosen 
maximum cruising power for all cases except that the 
effect of varying the cruising thrust has been deter- 
mined specifically for the jet airplane. Previous 
experience with conventionally powered airplanes 
has shown that the lowest operating costs are obtained 
with maximum cruising power usable with auto lean 
carburetor setting, since the reduction in maintenance, 
depreciation, and crew costs per mile obtained with the 
higher speed more than offset the higher fuel con- 
sumption. Maximum cruising power has been taken 
as two-thirds of the maximum continuous poweg for 
the DC-3 engine, which has a single-speed blower; 
two-thirds of maximum continuous power in low 
blower and about 69 per cent of that power in high 
blower for the DC-6; and about 85 per cent of maxi- 
mum continuous thrust for the jet transport. DC-6 
cruising ratings, 1,200 b.hp. per engine in low blower 
ratio up to 11,800, ft. and 1,100 b.hp. in high blower 
ratio, are those now in use by the air lines and are 
less than the maximum cruising powers that may be 
expected with the Pratt & Whitney CA15 engines 
after some experience is gained on these power plants. 
The jet engine cruising thrusts were selected from a 
study of existing information on gas turbine per- 
formance and operation. 

Cruising altitudes were selected to show the effect 
of cruising near the minimum and maximum altitudes 
practicable with the type of airplane being considered. 
The DC-3 has been studied at 1,000 and 10,000 ft. 
The DC-6 has been compared at 10,000 and 20,000 
ft., and the jet transport at 10,000 ft. and at its opti- 
mum cruising altitude, 35,000 to 40,000 ft. The jet 
transport achieves its greatest fuel economy at the 
highest possible altitude, as will be shown later. In 
order to permit sufficient airplane maneuverability, 
a practical operating ceiling was placed on the jet 
transport. The optimum altitude was therefore limited 
to the altitude at which the maximum rate of climb 
with cruising power is 350 ft. per min. The purpose 
of this limit is to avoid cruising too close to the mini- 
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mum drag speed, where it is difficult to maintain 
stabilized cruising speeds. Actually, a somewhat 
smaller rate of climb might be acceptable, but the 

cific value used iS not important to this study. 
The effect of this kind of limiting altitude is to make 
the cruising altitude a function of weight so that as the 
fuel is consumed, the airplane’s cruising altitude 
increases. 

The comparison of the three types of airplanes has 
been done at ranges considered typical for normal 
operation with these airplanes. The ranges used are 
400 miles for the DC-3, 1,200 miles for the DC-6, and 
750 miles for the jet transport. 


Descent 


Several types of descent have been studied with each 
airplane. In the case of the unsupercharged DC-3, 
these included descents with maximum cruising power 
at the maximum rate of descent permissible for pas- 
senger comfort and also at approximately the air speed 
for minimum sinking speed using whatever power is 
required to avoid exceeding this rate of descent. 

Another limitation that must be put on descents is 
a maximum value of indicated air speed. In the case 
of the DC-3, the low descent rate permissible prevented 
the air speed from exceeding the placard glide speed. 

For the pressurized DC-6, higher descent rates are 
possible. Therefore, in addition to studying the 400 
ft. per min. descents with maximum cruising power, 
descents were made with no power and with maximum 
cruising power at a.constant limiting indicated air 
speed. The limiting speed used was 300 m.p.h., 
indicated, since the placard glide speed of 360 m.p.h. 
was found to have excessive ride roughness, as well as 
to exceed the limiting Mach Number at altitudes above 
7,200 ft. Even the 300 m.p.h. indicated speed exceeds 
the Mach Number limit at the highest altitudes, as 
shown in Fig. 2 which gives the present speed limitations 
of the DC-6. A descent from 20,000 ft., limited by the 
300 m.p.h. speed and a cabin rate of descent of .400 
ft. per min. (1,000 ft. per min. for the airplane), was 
also studied, using whatever power was required up 
to maximum cruising power. The cases with constant 
300 m.p.h. indicated air speed exceed the allowable 
rate of descent but are included to show the effects on 
performance factors of such high descent rates. 

For the jet airplane, descents were computed at two 
constant indicated air speeds. One, 320 m.p.h., was 
chosen as the maximum speed at which good ride 
roughness qualities could be obtained. Actually, 
the speed might have to be reduced below this value 
at high altitude to avoid flight at an excessive Mach 
Number. To show the effect of reducing the descent 
speed, descents at 275 m.p.h., indicated, were studied. 


The descents have been computed with zero thrust 
and with one engine at maximum cruising thrust. The 
zero thrust case is of academic interest only, since the 
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turbines are used to supercharge the cabin and at 
least one engine must operate with appreciable thrust 
to supply the pressurization. 


Operating Costs 


Direct operating costs in terms of cents per two 
hundred pound-mile have been calculated by the 
standard method prescribed by the Air Transport 
Association (A.T.A.). In determining block time, 
the A.T.A. maneuvering time—0.12 hour for two- 
engined airplanes, 0.17 hour for four-engined airplanes— 
has been added to the climbing and cruising time. A 
10-m.p.h. headwind and reserve fuel adequate for 200 
miles plus 45 min. have been assumed in all range 
calculations. Fuel costs used were 13 cents per gal. 
of 91 octane gasoline for the DC-3, 16 cents per gal. 
of 100 octane gasoline for the DC-6, and 12 cents per 
gal. of kerosene for the jet airplane. 

The take-off weight limitation of the DC-3 deter- 
mines the pay load, since the nature of the structural 
design is such that the total of the weight empty, the 
fuel, and the pay load must not exceed this value. 
Because the fuel is placed in the inboard portions of 
the wing on the DC-3, the fuel does not contribute a 
relieving bending moment to the wing as it does in the 
DC-6 and the jet transport, where the fuel is largely 
in the outboard portion of the wing. The DC-3 is 
designed for a maximum take-off weight of 25,200 Ibs. 
and has an operating weight empty of 18,087 Ibs. 
The difference, 7,113 lbs., may consist of any combina- 
tion of fuel and pay load. Since full pay load and full 
fuel would considerably exceed this value, the allowable 
pay load is determined by the fuel requirement. For 
the DC-6 and the jet transport, the critical structural 
weight is the maximum permissible weight with zero 
fuel. All weight above thi zero fuel weight, which is 
70,000 Ibs. for the DC-6, for example, must be in fuel 
placed in the wing. These airplanes are designed 
so that they can carry almost full fuel above the zero 
Thus, the pay load is determined by the 
For 


fuel weight. 
zero fuel limit except for the longest ranges. 
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the ranges considered in this paper the full pay load 
can be carried. Another possible limit on pay load 
arises if the difference between the zero fuel weight and 
the landing weight is less than the required reserve fuel. 
With the reserves assumed in the studies discussed in 
this paper this condition does not arise. 


Ride Roughness 


Passenger comfort is affected not only by the rate 
of decrease of the cabin pressure altitude but also by 
ride roughness due to gusts. The ride roughness was 
specified by the magnitude of the accelerations due to 
sharp-edged gusts expressed in terms of the acceleration 
of gravity. The maximum indicated gust velocity of 
30 ft. per sec., used by the C,A.A. in establishing gust- 
load design conditions, was used as a basis of com- 
parison. The maximum acceptable ride roughness 
factor—i.e., the acceleration caused by a 30 ft. per 
sec. gust—was taken as 2.00, approximately the value 
experienced by a DC-3 in cruising conditions. 

The ride roughness factor is computed from the 
formulas 


An = KmU,V;/391 (W/S) 


and 
m = mo/(WV1 — M? cos* 0 + (19/xA.R.)] 
where 
An = ride roughness factor 
m = lift curve slope, dC,/da, per rad., corrected 
for aspect ratio, sweep, and compressi- 
bility effects 
mM, = incompressible section lift curve slope, 
dC,/da, per rad. 
U; = indicated velocity of the gust, m.p.h. 
V; = indicated airplane velocity, m.p.h. 
W = weight, Ibs. 
S = wing area, sq. ft. 
A.R. = wing aspect ratio 
M = = Mach Number 
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K = gust factor 
6 = angle of sweep of 25 per cent chord line 
of the wing 


The gust factor, K, adjusts the gust velocity to an 
effective sharp-edged gust velocity, KU;, by approxi- 
mately correcting for the lag in the development of 
lift on a wing when the angle of attack is changed 
suddenly and for the gradual build-up of acceleration 
of the airplane as the lift due to the gust develops, 
The use of the 30 ft. per sec. gust by the C.A.A. was 
a result of using a particular airplane for tests and a 
K of 1.0. Higher wing loadings cause larger accelera- 
tions than indicated by the equation, and thus K can 
be greater than 1.0. XK is given as a function of wing 
loading in Fig. 9 of reference 1. To correct for the 
effect of Mach Number on wing lift curve, the Glauert 
factor, 1/V 1 — M?, has been applied to the section 
life curve slope in computing the ride roughness factor, 
The Mach Number, M, is modified by the cosine of the 
angle of sweep, 6, to account for the reduced Mach 
Number effect on swept wings. One characteristic 
of sweptback wings which has been neglected but 
which is worthy of mention is the torsional flexibility 
of such wings. This would tend to alleviate some of 
the gust effect and would thus further reduce the ride 
roughness. 

The ride roughness factors for cruising and at the 
beginning and end of the descent are included in 
parentheses on the figures showing flight paths. The 
cruising values are taken at the end of cruising, since 
the most severe ride roughness occurs at the lightest 
weights. No ride roughness values are shown in 
climb, since the climbing speeds are relatively low and 
ride roughness would be less critical than in cruise 
and descent. 


RESULTS OF STUDY 


DC-3 Airplane 


Fig. 3 shows the effect of altitude on the operation 
of the DC-3 airplane. Comparing the cases using 
maximum cruising power for climb, cruise, and descent, 
the latter being at a rate of descent of 400 ft. per min. 
as are all the DC-3 descents, it is seen that the higher 
altitude increases the block speed by almost 7 m.p.h. 
and reduces the operating costs by 5.5 per cent. About 
one-third of the cost decrease is due to the increase in 
pay load possible because less fuel is required for the 
higher altitude, higher speed case. The ride roughness 
factor is reduced from 2.26 at 1,000 ft. to 2.11 at 
10,000 ft. 

Fig. 3 also includes a case for which rated power is 
used in climb. The block speed is increased by 5.5 
m.p.h. by this operation, but the costs are increased 
by about '/: per cent because the greater fuel con- 
sumed outweighs the cost advantage of the higher 
block speed. 
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Fig. 4 shows the effect of a descent made at a low 
indicated air speed, 120 m.p.h., with the least power 
required to keep the rate of descent from exceeding 
the 400 ft. per min. limit. This slow descent reduces 
the block speed by 14 m.p.h. and increases the cost 
about 3 per cent. Although the fuel required for the 
entire flight is reduced 4 per cent and the pay load 
increased correspondingly, the longer time required 
for the flight results in the higher costs, since all costs 
pased on hours of flying time are increased. The maxi- 
mum ride roughness, which occurs at sea level, is 
reduced from 2.54 in the maximum cruising power 
descent to 1.43 with the low air speed. 


From these cases it appears that the best form of 
operation from a purely performance standpoint is the 
high-altitude case using rated power for climb and 
maximum cruising power for cruising and descent. 
The use of rated power for normal climb would tend 
to shorten engine life and raise the sound level, how- 
ever, so that the recommendation might be modified 
to using the highest climb power compatible with 
satisfactory sound level and engine life. This type of 
operation would give the highest block speeds and 
very nearly the lowest direct operating costs. Under 
some conditions it might be necessary to reduce the 
descent speed to keep ride roughness within acceptable 
limits. 

It should be noted that with the DC-3 a constant 
take-off weight has been used for all cases and that 
the pay load is aways below the maximum volume 
capacity. Thus, any increase in fuel required calls 
for a decrease in pay load. If this airplane were 
permitted higher take-off weights so that the pay load 
could be held constant when fuel was added, operations 
that use more fuel would not show cost penalties due 
to decreased pay load. Then the rated power limit 
case, which has a lower pay load because more fuel is 
consumed in the climb, would show lower costs and 
would be the optimum case from both speed and cost 
standpoints. A generalization could then be made 
that the optimum operation is that using the highest 
power and altitude, provided the altitude does not 
exceed the engine critical altitude. 


DC-6 Airplane 


The various operations studied for the DC-6 airplane 
are plotted in Figs. 5 and 6. Fig. 5 shows the effect 
of cruising at 10,000 and 20,000 ft. and of using either 
maximum cruising or rated power in the climb. -All 
these cases use: maximum cruising power in cruise and 
descent and limit the airplane rate of descent to 400 
ft. per min. Fig. 6 illustrates the effects of four 


different descent procedures with all climbing done at 


rated power and all cruising at maximum cruising 
power. 

Fig. 5 shows that increasing the cruising altitude 
increases the block speed and decreases the cost for 
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both cruising power and rated power climb operation. 
For all these cases the pay load is constant, since the 
DC-6 is designed to carry full pay load with as high as 
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92 per cent of the maximum fuel capacity. For the 
1,200-mile range studied here, only about half of the 
normal capacity fuel load need be carried. Thus, 
changes in operating costs are due only to block speed 
and the cost of fuel consumed. One interesting result 
is that, if the climb is made at maximum cruising 
power, the gain in block speed resulting from flying 
at 20,000 ft. instead of at 10,000 ft. is only 5.6 m.p.h. 
or about 2 per cent, although the average cruising air 
speed is increased from 298 m.p.h. at 10,000 ft. to 
311 m.p.h. at 20,000 ft., a difference of 13 m.p.h. If 
the climbs are made at rated power, the block speed 
is increased 8.7 m.p.h. or 3.2 per cent. The reason 
for this effect is the greater time lost in climbing to the 
higher altitude with maximum cruising power. The 
shorter the time spent in climbing, the greater will 
be the net speed gain due to flying at the higher speeds 
possible at the higher altitude. Another reason for 
the relatively small difference shown due to altitude 
is that the maximum cruising power is reduced at 
20,000 ft. because operation in high blower ratio is 
required. The variation in speed with altitude is 
considerably greater with constant engine power, a 
condition existing when operating with the same blower 
ratio at altitudes below the full throttle altitude. 

Although there is a greater effect on block speed 
because of increasing the cruising altitude when the 
climb is made at rated power rather than at cruising 
power, the effect of cruising altitude on cost is inde- 
pendent of the climbing power. Costs are about 2.3 
per cent lower at 20,000 ft. than at 10,000 ft. in both 
cases. 

The comparison of operation at the same altitude 
with cruising or rated power in climb shows a block 
speed increase of 4.4 m.p.h. at 10,000 ft., due to utilizing 
the rated power climb, and a corresponding increase of 
7.5 m.p.h., when cruising at 20,000 ft. Thus the 
influence of climbing power on block speed is almost 
as great as that of cruising altitude. This is not so 
for the costs, however, since the cost reduction due to 
increasing the climbing power is only about '/2 per cent. 
The saving in block speed is almost canceled by the cost 
of additional gasoline. 


One advantage of the higher altitude is that the ride 
roughness factor is reduced from 1.68 at 10,000 ft. to 
1.50 at 20,000 ft. Both of these values are considerably 
less than the 2.11 experienced by the DC-3 at 10,000 ft. 

All the cases considered thus far have assumed 

descents with maximum cruising power at a rate of 
descent of 400 ft. per min. The maximum indicateg 
speed attained in these descents is nearly 300 m.p.h, 
and the ride roughness factor is always less than the 
DC-3 cruising value. With the pressurized DC4, 
however, greater rates of descent are permissible, 
Since the pressure differential between the atmosphere 
and the cabin is limited, the cabin altitude is a function 
of the airplane altitude. Fig. 7 illustrates this varia. 
tion, showing cabin altitude versus airplane altitude 
for the maximum permissible cabin pressure differential, 
From this relationship the minimum time of descent to 
sea level from any cruising altitude can be determined, 
assuming an average cabin rate of descent of 400 ft. 
per min. Fig. 8 converts this information into the 
maximum permissible rate of descent of the airplane 
from any cruising altitude with the cabin rate of 
descent limited to 400 ft. per min. The maximum 
airplane rate of descent from 20,000 ft. is 1,000 ft. 
per min., although a descent from a cruising altitude of 
12,000 ft. could be achieved at a rate of over 2,000 ft. 
per min. 

The effect of’ varying the type of descent has been 
studied using a 20,000-ft. cruising altitude as an 
example. To show the influence on performance, 
descents that exceed the 1,000 ft. per min. passenger- 
comfort limit have been included, since such rates of 
descent could be made usable by increasing the cabin- 
pressure differential and are possible even with the 
existing system when starting the descent from a lower 
cruising altitude. Another limit exceeded by the 300 
m.p.h. indicated air speed is the Mach Number of 0.54, 
which will permit a maximum indicated speed at 20,000 
ft. of only 278 m.p.h., as shown in Fig. 2. An indicated 
speed of 300 m.p.h. is permissible at 16,300 ft. 

Fig. 6 summarizes the results for four types of 
operation in the descent. All the cases assume a rated 
power climb to the cruising altitude, 20,000 ft., and 
cruising with maximum cruising power. Case 4, the 
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saine as that shown in Fig. 5, uses a rate of descent of 
400 ft. per min. with maximum cruising power. The 
indicated speed varies from 252 m.p.h. at 20,000 ft. 
to 291 m.p.h. at sea level. Case 5 keeps the same 
power but uses a constant indicated speed of 300 
m.p-h. For this case the rate of descent varies from 
1,685 ft. per min. at 20,000 ft. to 600 ft. per min. at 
sea level. Case 6 maintains the 300 m.p.h. indicated 
speed but descends with zero thrust, with the result 
that the rate of descent varies from 2,890 ft. per min. 
at 20,000 ft. to 2300 ft. per min. at sea level. Case 7 
combines the high indicated air-speed type of descent 
with the 1,000 ft. per min. limit. The descent is 
started with maximum cruising power and as high a 
speed as possible without exceeding the 1,000 ft. per 
min, rate of descent. The initial indicated speed is 
283 m.p.h. The indicated speed can be raised as the 
altitude decreases until 300 m.p.h. is reached at 9,000 
ft. Then the power is reduced from 1,100 to 875 
b.hp. per engine at sea level, the rate of descent being 
maintained at 1,000 ft. per min. and the indicated air 
speed at 300 m.p.h. 


The results show that the largest effect on operating 
cost does not exceed two-thirds of 1 per cent, and the 
largest effect on block speed is just under 3 m.p.h. or 
1 per cent. The optimum case is the cruising power 
descent at 300 m.p.h., indicated, but the detrimental 
effect on speed and cost of limiting the descent to 
1,000 ft. per min. is less than 1 per cent. Actually, 
the comparison between Cases 4 and 7 shows that a 
long power-on descent at a low rate of descent gives no 
loss in cost and a slight improvement in speed compared 
to the most rapid descent possible with a cabin rate of 
descent of 1,000 ft. per min. The ride roughness 
factors at the start of descent are about 2.20 for all 
cases except those with the 400 ft. per min. rate of 
descent, for which the ride roughness is about 1.70. 
However, for terrain clearance and traffic control 
reasons the low rate of descent operation may be unde- 
sirable. If low ride roughness and rapid descent are 
required, an operation with the indicated air speed 
restricted to less than about 250 m.p.h. and just enough 
power to limit the descent rate to 1,000 ft. per min. 
would be most desirable. The effect of descent speeds 
lower than the cruising speed is to decrease block speed 
and increase costs, as was shown in the low descent 
speed case for the DC-3. 


The results of the DC-6 study favor a rapid climb 
to the engine high blower critical altitude, cruising 
at maximum cruising power, and a power-on descent 
at a low rate of descent. The latter may be modified 
to reduce the time of descent by using the highest 
permissible rate of descent, the loss being a very small 
block speed reduction. If improved ride roughness 


is required, the power and air speed should be re- 
duced at the expense of a further reduction in block 


speed. 
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Jet-Powered Airplane 


The effects of cruising altitude and type of descent 
on a 750-mile operation are shown in Fig. 9. The 
altitudes compared are 10,000 ft. (Cases 1 and 2) and 
the altitude for the best cruising economy in terms of 
miles flown per pound of fuel (Cases 3 and 4). De- 
scents in this figure are at 320 m.p.h., indicated. Cases 
1 and 3 are with power off in the descent, and Cases 
2 and 4 are with one engine operating at maximum 
cruising thrust. There are two significant differences 
here as compared to the type of variation previously 
observed in the DC-6 study. One is that the cost per 
200 Ib.-mile is 37 per cent higher at 10,000 ft. than at 
the optimum altitude, whereas the effect of altitude 
was relatively small for conventionally powered 
airplanes. Another feature is that, whereas the cost 
reductions in the DC-3 and DC-6 cases were largely a 
direct reflection of an increase in block speed, the block 
speed is actually less at the high altitude for the jet 
transport. Thus, the cost reduction must result from 
a large decrease in fuel consumption, since all other 
costs, being dependent on block time, will be greater 
at the higher altitude. The ride roughness factor at 
the optimum cruising altitude for the weight at the end 
of cruising is about the same as for the DC-6 at 20,000 
ft. and is slightly lower at 10,000 ft. than for the DC-3 
at the same altitude. 


The effect of using power in the descent is to extend 
the descent over a longer time and distance. If the 
average air speed in the descent is higher than in 
cruising, the block speed is increased by using a low 
rate of descent. The average descent speed of the jet 
transport is lower than the cruising speed, so that a 
rapid descent with power off shows a higher block speed 
than does the power-on case. To permit average 
descent speeds as high as the cruising speed, the indi- 
cated speed at sea level would have to exceed 400 
m.p.h., thus imposing greater structural requirements 
and excessive ride roughness. The ride roughness 
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is about the same at 320 m.p.h., indicated, at high alti- 
tude as for the DC-3 in cruising. The ride roughness 
is higher at the optimum altitude and 320 m.p.h., 


indicated, for the zero thrust descent case than for the © 


partial power descent, because fuel for the descent need 
not be carried and the weight is therefore less. 


It is of considerable interest to note that the jet 
transport cruising ride roughness factor is lower than 
that of the DC-3 in spite of the greatly increased speeds. 
This is due to the increased wing loading, the low aspect 
ratio, and the sweepback of the wing, which reduces 
the increase in wing lift curve slope associated with 
high Mach Number. The influence of sweepback and 
aspect ratio on ride roughness is appreciable. The 
jet airplane considered in this study has a sweepback ~ 
angle at the wing 25 per cent chord line of 37.5° and 
an aspect ratio of 5.0. The ride roughness factor at 
the start of a 320 m.p.h. indicated air-speed descent 
from the optimum cruising altitude of Case 3 is 2.15 
with this sweepback and aspect ratio. Although the 
Mach Number is so high for this flight condition that 
sweepback would be necessary for satisfactory control, 


it is interesting to note that, if the sweepback were 
reduced to zero, the ride roughness would be increased 
to 3.18. If the aspect ratio of the jet airplane were 
equal to that of the DC-6 (9.45) and the sweepback 
angle were zero, the ride roughness for this condition 
would be 4.20. It therefore appears that ride rough- 
ness may be as important in selecting wing plan- 
forms for jet-powered airplanes as the usual charac- 
teristics such as compressibility drag divergence speed, 
induced drag, maximum lift coefficient, and wing 
weight. 

The effect of reducing the descent speed to 275 
m.p.h. is shown by Cases 5 and 6 in Fig. 10. These 
conditions are for maximum cruising power climb and 
cruise and use one engine at maximum cruising thrust 
in the descent. They compare to Cases 2 and 4, re. 
spectively, which differ only in the descent speed. As 
would be expected, reducing the descent speed de- 
creases the block speed by about 2 per cent for the 
10,000-ft. operation and twice that amount for the 
high-altitude cruising condition. Costs are increased 
by about 1 and 2 per cent, respectively. The initial 
rate of descent at the high cruising altitude, even at the 
reduced 275 m.p.h. descent speed, is 2,470 ft. per min. 
This high rate of descent would be made possible by 
the design of the cabin for the large pressure differential 
required for such a jet transport. From these results 
it appears that the optimum descent would be at the 
highest speeds permitted by structural and control 
requirements with as low thrust as possible. 

Fig. 11 shows the effect of using rated thrust rather 
than maximum cruising thrust in climb to be a 3 per 
cent increase in speed and a 21/2 per cent decrease in 
cost. The fuel used for the entire flight is actually 
less by several per cent when the climb is made at 
rated thrust. This differs from the conventional air- 
plane for which similar speed advantages are shown 
from using rated power in climb but for which most 
of the cost advantage derived from the speed improve- 
ments is lost in a higher fuel consumption for the trip. 
The reason for this is that the increase in the specific 
fuel consumption of the DC-6 engines at the high 
climbing power is not balanced by flying longer at 
20,000 ft. where only a 10 per cent increase in cruising 
miles per pound. of fuel is obtained with maximum 
cruising power in high blower as compared to the Fuel 
economy at 10,000 ft. in low blower. The turbojet- 
driven airplane obtains a large improvement in miles 
per pound at high altitudes and reduces its fuel con- 
sumption sufficiently by flying longer at the high 
altitude to overcome the eect of climbing with higher 
power. 

The variation of cost with power and speed at one 
altitude is of considerable interest. Since it was also 
desired to study the effects of throttling back at low 
altitude to reduce ride roughness, the 10,000-ft. alti- 
tude was selected. Cases 8 and 9 are identical to 
Case 6 except that cruising was done at the speed for 
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minimum drag-velocity ratio and at 10 per cent above 
that speed, respectively. The speed for minimum 
drag-velocity ratio, Vp,y, is the speed for maximum 
range if the specific fuel consumption is independent 
of thrust and speed as is roughly the case with turbojet- 
propelled airplanes. The results given in Fig. 12 show 
that the cost is nearly unchanged when the speed 


changes over a wide range. In this particular case 
the same cost is obtained when the average cruisiag 
air speed is decreased from the maximum cruising 
power value in Case 6 of 491 to 379 m.p.h. for 1.1Vpv 
operation. A further decrease to 345 m.p.h. increases 
the cost by less than 4 per cent. Thus the gain due 
to higher speeds which reduces all the costs dependent 
upon flight time is almost exactly compensated by 
the higher fuel costs at the maximum cruising power 
speed. For the 10,000-ft. high power case the fuel cost 
is about 69 per cent of the total cost, whereas for an 
airplane like the DC-6 it would be approximately 
one-third of the total. At high altitudes the per- 
centage of the total cost chargeable to the fuel is re- 
duced to 56 per cent, so that gains obtained by increas- 
ing power and speed would show some reduction in 
operating costs. 

The greatly increased fuel consumption of the jet 
transport is illustrated by Fig. 13, which shows miles 
flown per pound of fuel for the DC-3, the DC-6, and 
the jet transport. The DC-6 carries over three times 
as much pay load as the DC-3, so that the increase in 
fuel per mile for the DC-6 is less than the increase in 
useful work performed. In this study the specific 
fuel consumption of the jet airplane was reduced 
10 per cent from that of existing engines to allow for 
future development. Even then the jet transport 
consumes two to four times as much fuel as the DC-6, 
although it carries a smaller pay load. The signifi- 
cance of the fuel consumption and the cost of the fuel 
on the operating efficiency of the jet airplane cannot be 
overemphasized. 

The increase in fuel required at low altitude has even 
more serious implications than indicated by the in- 
creased costs. The total fuel burned when cruising at 
the high altitude, as in Case 7, is 17,250 Ibs., and the 
fuel carried including reserve is 24,200 lbs. Using 
maximum cruise power to cruise at 10,000 ft., as in 
Case 6, the fuel actually burned is 27,550 Ibs. and the 
total fuel carried is 45,300 Ibs. The particular jet 
transport studied had a design take-off weight of -ap- 
proximately 76,000 Ibs., the gross weight obtained with 
a full pay load of 12,000 Ibs. and capacity fuel of 24,500 
lbs. Thus, the 10,000-ft. cruising case requires almost 
twice the fuel that can be placed in the airplane struc- 
ture. Although some increase in fuel capacity can be 


obtained by the installation of tip tanks, the large 
additional amount of required fuel is clearly excessive. 
Even if space could be found for over 45,000 Ibs. of fuel, 
the effects on take-off distance are rather devastating. 
The take-off weight is increased from 76,000 to 96,800 
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Ibs. and the corresponding C.A.A. minimum take-off 
runway lengths from about 5,100 to 9,500 ft. Reducing 
the speed at 10,000 ft. to the 1.1Vp,;y speed would re- 
duce the total fuel carried to 37,500 Ibs., but even 
then the take-off weight, 89,000 Ibs., and the take-off 
runway length, 7,600 ft., are considerably above the 
design values. If the airplane were actually to be 
designed for these high weights, some increase in 
structural weight would also be necessary. Thus, 
because of lack of fuel capacity, excessive take-off field 
lengths, and excessive costs, operation at low altitudes 
would be prohibitive. The jet transport is literally 
forced to operate at the highest possible altitudes if 
reasonable ranges are to be obtained. 

The detrimental effect of low-altitude operation 
must be considered in specifying reserve fuel. In all 
the jet transport studies, the reserve adequate for 200 
miles plus 45 min. of flight was assumed to be burned 
at the cruising altitude and power. For the con- 
ventionally powered airplanes the reserve has been 
assumed to be used at 10,000 ft. or the cruising altitude, 
whichever was the lower, and at a speed 125 per cent 
above the maximum range speed as specified by the 
A.T.A. method. The effect of varying the reserve 
fuel on the DC-6 is to slightly vary the average flying 
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weight, since the reserve fuel is only 3.4 per cent of 
the zero fuel weight and full pay load could be carried 
in all cases. The DC-3 pay load would be affected 
by the reserve fuel that must be carried, since the 
take-off weight is always limiting for the cases studied, 
“but here, too, the effect is small, since the reserve fuel 
is only 3.7 per cent of the zero fuel weight. Further- 
more, the altitude at which the reserve fuel is consumed 
has a relatively small effect on the quantity of reserve 
fuel required for conventionally powered airplanes. 
The jet transport reserve fuel is a large quantity, 
however, being 13 per cent of the zero fuel weight when 
the reserve is burned at high altitude at maximum 
cruising power and 34 per cent of the zero fuel weight 
when the reserve fuel is burned at 10,000 ft. with maxi- 
mum cruising power Even when the speed at which 
the reserve is used at 10,000 ft. is reduced to 10 per cent 
above the speed for minimum drag-velocity ratio, the 
reserve fuel is 24.5 per cent of the zero fuel weight. 
Thus, decreasing the reserve fuel altitude from the 
optimum cruising altitude to 10,000 ft. will increase 
the fuel that must be carried by 11 to 21 per cent of 
the weight of the airplane with zero fuel, depending 
upon whether a low or high cruising thrust is used at 
the 10,000-ft. altitude. This corresponds to a weight 
change on this particular airplane of 5,800 to 11,000 
Ibs., which would cause appreciable increases in the 
required take-off and landing field lengths and an 
increase in fuel capacity required of 24 to 45 per cent— 
an increase that would be extremely difficult to realize. 
It is therefore clear that the jet transport must be so 
operated that as much as possible of its reserve fuel is 
consumed at the optimum altitude. This fact intro- 
duces an entirely new type of airways traffic control 
problem, because such an airplane would not be able 
to circle an airport at altitudes of 5,000 to 10,000 ft. 
for long periods. 


Although detailed attention has not been given to 
short-range operation with the jet transport, it can be 
seen from Fig. 9 that for ranges less than about 400 
miles the most efficient operation of the airplane would 
tend to consist entirely of a climb and a descent. This 
would introduce additional traffic control problems, 
since such a flight plan would not even have such a 
thing as a cruising altitude. Considerable study would 
be required to determine a satisfactory method of 
safely regulating airways populated by airplanes that 
are almost always climbing and descending and seldom 
flying level. The existence of optimum altitudes that 
greatly improve airplane operating efficiency will also 
tend to crowd certain altitude regions. 


The reasons for the greater effect of altitude on the 
jet transport than on the propeller-driven airplanes 
are to be found in the nature ofthe turbojet-powered 
airplane. There are three types of changes that affect 
the fuel economy of any airplane when the altitude is 
(1) Consider the angle of attack, and thus 


altered. 
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the indicated speed, to be held constant as the altitude 
is increased. The level flight speed then increases as 
the reciprocal of the square root of the air-density 
ratio. The drag is unchanged so that the thrust 
horsepower required rises with the velocity. Assuming 
that for the propeller-driven airplane the propeller 
efficiency and the specific fuel consumption (pounds 
of fuel per hour per horsepower) are constant, the miles 
flown per pound of fuel is unchanged. For a constant 
thrust unit, the miles flown per pound of fuel increase 
with the velocity, assuming the specific fuel con- 
sumption (pounds of fuel per hour per pound of thrust) 
to be constant. The essential difference is that the 
fuel used per hour by the propeller-driven airplane is 
dependent approximately on the product of drag times 
velocity, while the jet airplane fuel consumption per 
hour depends only on drag. Thus, to improve fuel 
consumption per mile, the drag of the propeller-driven 
airplane must be reduced, while the jet airplane gains 
from either a drag reduction or a speed increase. (2) 
When maximum cruising power operation is being con- 
sidered, the power is constant for the conventional 
engine as the altitude increases up to the critical 
altitude, while the thrust of the jet engine decreases. 
In either case the airplane cannot be flown at the higher 
speed necessary to maintain constant angle of attack. 
The angle of attack therefore increases, decreasing the 
drag. This benefits the propeller-driven airplane 
directly, and, while the speed decrease counteracts 
some of the gain in miles per pound of fuel for the 
jet airplane, the increase in the angle of attack brings 
the airplane closer to the maximum range attitude for 
constant thrust airplanes—i.e., the attitude for mini- 
mum drag-velocity ratio. (3) The third factor affect- 
ing range when the altitude is increased is the variation 
in specific fuel consumption. 


Fig. 14 shows the percentage change with altitude 
of the specific fuel consumption with maximum cruising 
power for the three airplanes considered in this paper. 
The DC-3 shows a constant specific fuel consumption 
because it operates in one blower ratio at maximum 
cruising power with a constant engine speed at all 
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altitudes. The DC-6 has essentially the same charac- 
teristics, except that the specific fuel consumption 
decreases as the power drops at altitudes above the low 
blower critical altitude and increases when the blower 
ratio is changed. These data assume’ perfect. car- 
buretor altitude compensation. The jet engine shows 
a large decrease in specific fuel consumption as the 
altitude is increased. This increase in fuel economy 
is due to the decrease in intake air temperature and is 
characteristic of the gas turbine. 

Fig. 15 shows the percentage change with altitude of 
the miles flown per pound of fuel. The DC-3 and the 
DC-6 show the gain in fuel economy due to flying at 
an angle of attack with less drag, the higher air speed 
in itself, and the specific fuel consumption having no 
effect. The DC-6 shows a reduction in the miles 
flown per pound of fuel which occurs when the blower 
ratio is changed. The jet transport gains from all 
three of the factors described above, so that the increase 
in miles per pound is extremely large. The increase in 
speed, the improved angle of attack, and the improved 
specific fuel consumption combine to make the miles 
flown per pound of fuel more than 100 per cent greater 
at 35,000 ft. than at 10,000 ft. 
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CONCLUSION 


This study of the effects of various operational param- 
eters of transport aircraft, such as engine power, 
cruising altitude, and rate of descent, has shown that, 
for both conventionally powered and turbojet-powered 
airplanes, the best: economy and highest block speeds 
are obtained by climbing rapidly, cruising at high alti- 
tude, and using as much power as practicable. Be- 
cause the descent with a DC-3 or DC-6 type of airplane 
can be flown at an average true speed above the cruis- 
ing speed, it is desirable for best speed and economy to 
extend the descent time by using cruising power. The 
jet airplane operates at such high altitudes and speeds 
that average descent speeds as high as the cruising 
speeds cannot be attained with satisfactory ride rough- 
ness. The fuel consumption is also excessive at low 
altitudes. For best jet airplane speed and economy it 
is then desirable to descend as rapidly as possible with 
the minimum possible thrust. Cabin pressurization 
and ride roughness may restrict rate of descent and 
speeds, respectively, even if structural and control speed 
limits are not exceeded. 

While the advantages to be gained froin operating 
the airplane in the optimum manner are relatively small 
for conventionally powered aircraft, the jet airplane 
must be operated at high altitude to achieve reasonable 
costs and ranges. These quantities are so severely 
affected by cruising at low altitudes that the jet trans- 
port becomes relatively inefficient, if not impractical, 
for such operation. Because of the extremely high 
altitudes that are required for efficient operation, even 
for short ranges where the flight path reduces to only a 
climb and descent, the problem of developing a satis- 
factory system of airways traffic control for jet trans- 
ports is a serious one. 
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The Effects of Compressibility on the 
Two-Dimensional Subsonic Wind-Tunnel 
Constriction Correction 


H. T. EPSTEIN* anp L. U. ALBERS? 
Unwersity of Michigan 


SUMMARY 


Relaxation methods were applied to the determination of the 
effects of compressibility on two-dimensional wind-tunnel iuter- 
ference. Calculations were made for a 10-per cent-thick sym- 
metrical section, both isolated and symmetrically placed between 
two walls, for ratios of chord to tunnel height c/h of 1 and 2. 
Flows at Mach Numbers of both 0 (incompressible flow) and 
0.5 were obtained. 

The incompressible velocities agreed well with those obtained 
by conformal mapping methods and the compressible velocities for 
the free-stream case were in satisfactory agreement with those 
given by Kaplan.‘ 

Evaluation of the wall interference showed that the effect of 
compressibility averaged linearly along the chord was about the 
same as given by the linear-perturbation theories of Goldstein 
and Young, and Allen and Vincenti. The variations along the 
chord were so great, however, that the use of these theories for 
correcting local quantities is problematical. 


INTRODUCTION 


4 Beer OF LARGE MODELS at subsonic compressibility 
speeds in wind tunnels have necessitated theoretical 
investigation of the effects of compressibility on wind- 
tunnel corrections. An important correction is that 
of constriction, which gives the effect of the finite cross- 
sectional area of the model in constricting and acceler- 
ating the flow. 

Investigations of the constriction effect have been 
made by Goldstein and Young! and by Allen and Vin- 
centi.? The results of both investigations show that 
incompressible velocity increments due to the tunnel 
need only be multiplied by (1/ V1 — M,°)* to obtain 
the compressible velocity increments, where M) is the 
Mach Number in undisturbed stream. Both theories, 
however, assume small thickness-ratio airfoils and small 
values of the ratio of chord to tunnel height o. In 
addition, the result is an ‘‘average” correction in the 
sense that (in the Goldstein-Young analysis) the air- 
foil is replaced by an equivalent doublet and the re- 
sultant correction applies as an overall or average cor- 
rection. Thus, the applicability of these formulas to 
large models (large o values), to cases of large thickness 
ratios, and to local quantities on the airfoil had to be 
investigated. The evaluation of the constriction cor- 
rection for a symmetrical flow is given in this paper, 

Received September 15, 1947, 


* Department of Physics. 
+ Department of Mathematics. 


144 


which considers a thickness ratio of 10 per cent and 
o values up to 2. 

Among the theoretical methods of calculating in- 
-compressible and compressible flow patterns, relaxation 
methods appear promising in view of the relative ease 
with which boundary conditions may be applied. The 
principle of the relaxation solution of partial differen- 
tial equations is to satisfy in a systematic manner 
the finite difference form of the equations at a pre- 
scribed number of net points of the field.* 

Velocity calculations were made for Kaplan’s 10 per 
cent thick symmetrical section,‘ both isolated and for 
o ratios of 1 and 2 and for stream Mach Numbers of 
0 and 0.5. From these results, the compressible and 
incompressible velocity increments due to the tunnel 
were obtained. The values are compared with those 
predicted by Goldstein and Young', and Allen and 
Vincenti.? (See Figs. 1 to 4.) 

The Kaplan section was chosen for the following 
reasons: 

(1) Both the compressible and incompressible free- 
stream flows had already been calculated by other 
writers, so that comparison of results was possible. 

(2) It is easily possible to verify the theoretical 
predictions in any of the small tunnels now built in this 
country (see section on “Comparison with Experi- 
ment’’). 


SYMBOLS 


cross-sectional area of flow 

velocity of sound 

ratio of airfoil chord to tunnel height (c) 

Mach Number (V/a) 

exponent of 8, defined by K = 1/,” 
m + m2 + 93 + — 40 

arbitrary constant 

velocity 

V1 — Me 

Vi Vr 

incompressible stream function 
angle surface makes with x-axis 
incompressible velocity potential 
density 
compressible stream function 


Subscripts 


0, 1, 2,3, 4 = points of net 
= compressible 
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i = incompressible 

= isolated 

0 = undisturbed stream 

ST = stagnation condition 

x, ¥ = partial derivatives in the x and y directions 


RELAXATION METHODS OF CALCULATION 


The procedures given by Emmons for performing 
relaxation solutions for compressible and incompressible 
flows are used in this paper. A detailed discussion of 
the relaxation methods is given in the appendix; only 
an outline of the compressible procedure is presented 
at this point. 

In this treatment of a compressible flow by relaxation 
methods, it is necessary to know the incompressible 
flow past the same body. The physical plane is trans- 
formed into a plane whose coordinates are the incom- 
pressible velocity potential and stream function. The 
difference equation for the compressible stream func- 
tion in terms of these coordinates is then systematically 
solved. The Mach Number at each point can be ob- 
tained from the stream function. 


COMPARISON METHODS OF OBTAINING VELOCITIES 


Incompressible Velocities 


The incompressible velocity distributions for the 
various tunnel configurations were supplied by Perl and 
Moses. The method used consisted basically in con- 
formally mapping the section into a lattice of straight 
lines. The mapping was accomplished numerically by 
an iteration method. 


Compressible Velocities 


Free-Stream Flow.—Kaplan‘ gave the free-stream re- 
sults by a successive approximation method similar 
to that of the Poggi method, except that the compres- 


sible stream function was expanded in powers of air- 
foil thickness ratio instead of Mach Number. 


Wind-Tunnel Corrections 


Using the Prandtl-Glauert rule as a basis, Goldstein 
and Young! considered the corrections to the velocity 
given by the simple Glauert image theory of tunnel 
corrections. By taking into account the linear trans- 
formations of coordinates and velocities called for in the 
Prandtl-Glauert analysis, they gave formulas for the 
effect of compressibility on wind-tunnel corrections in 
terms of powers of 8 where 8 = Vi- M,?. For the 
constriction correction, they gave 


K = (AV/V))-/(AV/ Vo); = 1/83 
RESULTS AND DISCUSSION 


Velocity Distributions 


The theoretical velocity distributions obtained with 
the Kaplan airfoil are given in Fig. 1. The correspond- 
ing results given by the comparison methods of refer- 
ences 4, 5, and 7 outlined in the preceding section are 
also indicated. 

The derived wind-tunnel-constriction corrections 
AV/Vo, the velocity increment due to the tunnel 
divided by undisturbed stream velocity, are presented 
in Fig. 2. The effect of compressibility on the correc- 
tion is shown in Fig. 3. The value from references 1 and 
2 is also indicated in Fig. 3. It should be noted that ¢ = 
0.25 has been considered to be equivalent to ¢ = 0 in 
Fig. 1. The agreement with the results of references 
4 and 7 is considered sufficient evidence for the correct- 
ness of the assumption. 

In general, the agreement of the relaxation methods 
with the other methods discussed is good. In Fig. la, 
the error near the leading edge in the incompressible 
free-stream flow is ascribed to the fact that only ten 
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50 
points on the airfoil were used in the relaxation, which io YA 
means that the effect of the leading edge of the cusp Sas B 
was given only by the leading-edge point itself. Thus, S 
for example, it would be difficult, to this accuracy, to a J 
distinguish between a cusp and any other airfoil with ral as 
a pointed leading edge. Because the percentage differ- 4 4 
ence in ordinates would be large, the results should be a 35 
correspondingly different. Inasmuch as greater ac- 
curacy was reyjuired for the problems with the airfoil in 
a tunnel, the accuracy of the results was improved by z / 
subdividing the first rough net into a finer net of points S / 
near the airfoil. a 25 
The compressible relaxation results in Fig. la are 8 7 ane ta 
uniformly slightly lower than the results of reference 4. 6 20 OD OF! 
This effect is ascribed to the fact that in the determina- & M —- AVERAGE OF RELAX- 
tion of the incompressible potential lines in the relaxa- 
tion method / (Oy/dy) dx was evaluated instead of 2 V4 —-- REFERENCES | AND 2 
JS ~/0 n) ds, where n refers to the direction of the Fa T 
normal to the surface of the airfoil and s to the direc- 40 BO BS 


tion of the tangent to the surface. The difference be- 
tween the two integrals was not observable in the 
graphical integrations performed, but the difference ob- 
viously exists. In fact, the difference is a factor 
of 


(ds/dx)* = 1 + tan? 6 


If the value of tan @ may be considered approximately 
equal to the thickness ratio, the factor has the value 
1.01. This 1 per cent adequately accounts for the dif- 
ference between the results of reference 4 and those of 
the present report. The factor is also omitted in the 
other relaxation results presented but does not affect 
the comparative value of the results. 


PERCENT HALF —CHORD 


Fic. 3. Effect of compressibility on constriction correction. 
Fic. 3a (top). K = (AV/Vo)comp./(AV/Vo)incomp.. Fic. 3b 
(bottom). n, defined by (1/+/1 — M,?)"= K. 


Wind-Tunnel Corrections 


Figs. 3a and 3b give K and 2, respectively, as func- 
tions of distance along the chord, where n is defined 
by 


K = 1/8" 


Because the results near the leading edge are uncer- 
tain, the curves have been faired in to make them 
continuously decreasing. This fairing is justified by the 
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GCy _N 


— 


Fic. 4. Net configurations for. 10 per cent thick section. 
Fic. 4a (top). Ten per cent thick section symmetrically placed 
in a o = 2. Fic. 4b (bottom). Relaxation net used for 
¢=2,M = 0. 


fact that a one-dimensional treatment shows that K 
and m decrease continuously toward the leading edge, 
and it is expected that one-dimensional theory gives 
the correct qualitative behavior. The reason for fairing 
is that an average correction was desired for comparison 
with the value m = 3 in references 1 and 2. For this 
purpose, the exact curve is not needed because the con- 
tribution of the faired part to the total area under the 
curve is not large. 

In Fig. 3b, the relaxation m and the value » = 3 from 
references 1 and 2 are indicated together with the re- 
laxation » averaged along the chord. It is seen that: 
(1) m varies widely along the chord; (2) the value » = 
3 is close to the linear average m of the relaxation re- 
sults; and (3) the value of m may be much higher and 
much lower than 3. Inasmuch as m is much higher 
than 3 at the point of maximum velocity, the effect of 
compressibility is to bring on phenomena (such as 
choking) very quickly as Mo is raised and thus limits 
the useful Mach Number range of a tunnel. 

Vincenti and Graham* examined the constriction 
correction for airfoils in closed circular tunnels and 
found the same m = 3 correction. Examination of cer- 
tain experimental data cited in reference 8 showed that 
the theory gives good results up to ¢ = 0.625. Thus, 
the present report may be considered to extend the 
range of validity of the formulas to o = 2 if the linear 
average of m along the chord can be identified with the 
average m given in references | and 2. 


Comparison with Experiment 


Owing to the symmetry of the Kaplan airfoil, it is 
possible to obtain the experimental flow in either a 
flexible wall tunnel or in a regular tunnel wherein one 


side has been modified to the correct shape. It can be 
seen from Fig. 4a that, since ADM is a streamline, it is 
possible to replace it by a solid wall and thus simulate 
the flow past the airfoil by the flow in a channel with 
walls ADM and BN, wherein the airfoil has become 
a bump in the wall. 

If this is done, however, it is necessary to correct 
for the boundary layer on the walls, and it is suggested 
that this can be done by the method of von Karman 
as presented by Puckett. Puckett presents formulas 
expressing the displacement boundary-layer thickness 
in terms of the actual boundary-layer thickness and 
the local Mach Number. The Mach Number can be 
obtained from pressure measurements, and the actual 
boundary-layer thickness measured from schlieren pho- 
tographs, so that the displacement thickness can be 
calculated. Then the Mach Number for zero boundary 
layer can be obtained from the ratio of the actual to the 
“effective” cross-sectional area. The effective cross- 
sectional area is the actual area minus the displacement 
thicknesses on the tunnel walls. _ 


CONCLUSIONS 


(1) The effect of compressibility on the two-dimen- , 
sional wind-tunnel-constriction correction varies con- 
siderably along the chord. The Goldstein-Young and 
Allen-Vincenti result of constant compressibility effect 
is approximately the same as the linear average effect 
of the calculated result. 

(2) The relaxation methods of obtaining two-dimen- 
sional incompressible and compressible subsonic flows 
give satisfactory results, inasmuch as they agree closely 
with those of other theoretical methods. 


Appendix 


THE RELAXATION PROCEDURE 


The general tunnel configuration is indicated in Fig. 
4a. Asa result of the two-way symmetry of the prob- 
lem, it is possible to consider only the quadrant ABCD, 
as shown in Fig. 4b. All numerical values given are 
taken from the solutions for ¢ = 2. 


Incompressible Flow 


The incompressible stream function » is taken to be 0 
on the streamline AED in Fig. 4b. For convenience, 
the top wall BC is called the 100,000 streamline. Inas- 
much as uniform flow is assumed to prevail at AB up- 
stream of the airfoil, the values of n at points a), de, etc., 
are simply 20,000, 40,000, etc. Along DC, the first 
approximation is that 7 increases linearly. Thus, the 
values at the net points are taken to be 25,000, 50,000, 
and 75,000, as indicated. Because there are no net 
points at the airfoil surface, it is necessary to find the 7 
values at the points nearest the airfoil and to kéep these 
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“boundary” values fixed instead of the airfoil-surface 
values, which means fixing the value of » at a point 
like b. The method of fixing the value will be 
explained shortly. All » values on the net-point 
boundaries are then held constant throughout the re- 
laxation. 

Initial values of » at all remaining points inside the 
net boundary are guessed by sketching in the stream- 
lines freehand and interpolating to the net points be- 
tween the streamlines. 

Incompressible flows are given by the Laplace equa- 
tion; 


Nez + = 0 
or its finite difference form :* 
m + n2 + 13 + m — 4m = 0 (1) 


where the points 0, 1, 2, 3, and 4 are as shown in Fig. 
5a. If the values of 7 at these points are incorrect, the 
left-hand side of Eq. (1) will not add up to 0 but will 
give a residual Q. ‘The purpose of the relaxation 
method is to reduce (or relax) the residual until it is 
sufficiently close to 0 for the difference to be negligible 
in comparison with the difference in » between neigh- 
boring net points. The residual Q at each point is thus 
calculated and recorded beside the 7 value at the point. 
(See Fig. 5b, which is the portion UVWZ indicated in 
Fig. 4b. Then, if the » value at a point such as X is 
changed from 53,422 to 53,423, the residual at that 
point is reduced from 3 to —1, whereas the surrounding 
residuals are changed by only the one unit by which mo 
was changed. Thuis, at each point the residual may be 
reduced to almost 0. It is wise to relax the point in the 
field with highest residual rather than to clear up a 
whole area, inasmuch as the effects of changes in one 
area spread out and may completely nullify the work 
of relaxing any other area. 

The symmetry of the problem permits the line CD to 
be effectively a boundary although the 7» values along 
it are not known. (See Fig. 4a.) The 7 values along 
GH are the same as along JK. If the 7 values on these 
lines are kept the same, the effect of the symmetry on 
the n values along CD is obtained. 

Reference is again made to Fig. 4b; the net point 
b is set by interpolation between the surface S and the 
point C, whose value was estimated. As the relaxation 
proceeds, the value of » at c changes. The change is 
eventually so great as to make the value of b incorrect; 
b is therefore recalculated by linear interpolation using 
the new value of c. This redetermination ordinarily 


changes the residual at c by such an amount that the 
frequent recalculation of 6 is recommended. After 
some experience, however, it is possible to judge the 
change in b so closely that only three of four changes 
are needed before the final values for the net points are 
determined. Velocities are then found from the for- 
mula 
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-1155,792 -254,596 3)/53,422 . 66 -1/51,518 
- 3,423 


-133,763 31,863 29988 -228289 -2 26,898 
(w) (2) 


Fic. 5. Portion of net for 10 per cent thick section; o = 2, 
M = 0. Fic. 5a (top). Basic net point location for difference 
equation. Fic. 5b (bottom). Portion of M = 0,0 = 2 net with 
one indicated relaxation. 
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where V, is the velocity at the point 0 and 6 is the 
stream-function difference between two net points in the 
undisturbed stream. 


Compressible Flow 


The assumption that the values of the stream func- 
tion near the surface may be set by linear interpolation 
is shown by Emmons* to be inadequate when com- 
pressible flows are considered. As a result, it is not 
ordinarily possible merely to use the compressible 
equation of flow in finite difference form and to proceed 
as outlined for the incompressible case. Emmons 
recommended, and the authors used, a scheme that 
consists in mapping the physical plane into the plane 
where the coordinates are the incompressible velocity 
potential and stream function. Since the compressible 
and incompressible potential and streamline patterns 
are reasonably similar, the required relaxation is rather 
small. In this plane the airfoil becomes a segment of 
one of the streamlines. The equation of irrotational 
flow is then expressed in terms of the new variables and 
in finite difference form is given by Emmons (with his 
w = 0)as 


Vi + vo tvs + — Sho — 


)| (ve. 57), | 


Fic. 6 


The b 
Numi 
lem c 
to be 
OF 
in 
surfac 
found 
arran 
Alo 
the st 
value 
wall i 
pressi 
est m 
the 7 
sumit 
streal 
sonic 
No 
much 
the m 
illusts 


4 
3 
(u) (v) 
-1]77305 -2]76,741 75,832 
-| 
is 
The 
power! 
Use 


o = 2, 
ference 
et with 


func- 
lation 
com- 
Ss not 
ssible 
mons 
that 
ocity 
sible 
terns 
ither 
it of 
onal 
and 
1 his 


(3) 


SUBSONIC WIND-TUNNEL CONSTRICTION CORRECTION 149 


Los =— 
Pst 
.630-M/ 638 641 65! 656 
-192 197 188 - 
1.182 1201 (5 1.226 
17276 17276 17276 (7276 
648 655 660 
-202  |-.206 
1207 1.234 
12957 12957 12957 
652 661 
-205  _|-.210 
1.220 |! 237. 
66! 669 675 
1.183 L213 1238 1.269 
4319 4319 4319 4319 4319 
Fic. 6. Portion of initial net for 10 per cent thick section; o = 


= VO. 


It is further shown that the nondimensional mass flow 
is 


pV _R 
= [Ws — + — (4) 
Pst 26 ‘ 
The ratio R/é is ordinarily taken to be some integral 


power of 10. 
Use is also made of the well-known relations 


oy = (5) 


Pstasr 


6 
log. OST log, [(: + (6) 


+1)] 


The boundary values can be set by specifying the Mach 
Number of the undisturbed stream, which in the prob- 
lem considered was 0.5. Eq. (5) then gives pV/psrasr 
to be 0.4319. Eq. (4) with R/S = 10- gives (Yi — 
v2) or with on the surface streamline, since = 
%, in uniform flow and since Y is taken to be 0 at the 
surface; thus, ¥2 = 4,319. The term log, (p/pST) is 
found from Eq. (6). The necessary quantities are 
arranged as shown in Fig. 6. 

Along the potential line far upstream of the airfoil, 
the stream function again increases linearly, and the y 
values are 8,638, 12,957, etc. Inasmuch as the top 
wall is a streamline, y is there 21,595. 

It is then necessary to estimate values of the com- 
pressible stream function throughout the net; the easi- 
est means is to take a constant value of y along each of 
the 7 lines. Keeping y constant is equivalent to as- 
suming that the incompressible and compressible 
stream lines coincide to a first approximation in sub- 
sonic flow. 


Now, consider any point such as X (Fig. 6). Inas- 


much as all the quantities of Eq. (4) are known except 
the mass flow, p V/psragr may be calculated to be, in this 
illustration, 0.527. Eq. (5) then gives M at that point 


to be 0.652. From Eq. (6), log, (p/pST) is found to 
be —0.205. When all these quantities are found for 
each net point, the residual can be determined at each 
point from Eq. (3). 

For nearly parallel flows, as for airfoils of small 
thickness ratios, the last term in Eq. (3) is negligible, 
and the next to the last term varies so slowly that, for 
a while at least, only the first five terms in the relaxa- 
tion need be considered; but these terms are precisely 
the incompressible ones. The relaxation therefore pro- 
ceeds as in the incompressible case until the quantities 
have so changed that the error in neglecting the last 
two terms is no longer negligible. It is then necessary 
to go through the previous procedure to obtain the 
quantities required for finding the residuals, this time 
using the newly found values of y. Again the relaxa- 
tion continues as on an incompressible problem; the 
correct residuals are occasionally calculated until the 
y values have converged. In the cases considered, this 
process converges after two or three recalculations of 
the residuals. 


In this fashion the Mach Number at each. point is 
found and from it the pressure and velocity can be de- 
termined by well-known relations. 


General Remarks on Relaxation Solutions 


In problems in which symmetry does not permit the 
reduction of the field, a further variable is present which 
is the circulation around the airfoil. The Kutta- 
Joukowsky condition, that the airfoil streamline leaves 
the airfoil at the trailing edge, is then used. The value 
of the airfoil stream function is estimated; the net is 
relaxed; and, if the streamline does not leave the sur- 
face at the trailing edge, the estimate is altered and the 
rather slight amount of relaxation required is per- 
formed. This procedure is also one that converges 
rapidly. Account can be taken of boundary-layer and 
separation effects simply by requiring that the airfoil 
stream function leave the surface at some point up- 
stream of the trailing edge. 

Velocities may be determined by proceeding in one 
of two ways. The first is to make the net finer and 
finer near the airfoil until the points are so close to the 
surface that the error is negligible in taking the velocity 
off the surface to be equal to that on the surface. The 
alternative is to plot the velocity along some vertical 
line and extrapolate from the net point nearest the sur- 
face to that value which would obtain at the surface 
itself. Both methods have been used in this investiga- 
tion with equivalent results. 
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Letter to 


Dear Sir: 

The article “Aerodynamic Performance of Delta Wings at 

Supersonic Speeds’’ by A. E. Puckett and H. J. Stewart (Jour- 

* NAL OF THE AERONAUTICAL SCIENCES, Vol. 14, No. 10, p. 567, 
October, 1947) contains, on p. 577, a method of determination of 
a suction force at the sharp leading edge of an airfoil when the 
edge sweep is such that the relative velocity component, normal 
to the edge, is subsonic (“subsonic edge”’). 

The writer feels that this force does not exist and that the 
method of calculation presented in the article does not apply to 
compressible fluids. 

In the case of incompressible fluids, the suction force T is (cf. 
Durand, W. F., Aerodynamic Theory, Vol. II, p. 52, referenced 
hereafter as A. T.) : 

T = xpC? (1) 
with 
u= +C/V/x (2) 
where u is the velocity perturbation normal to the edge and «x is 
the distance to the edge. This expression for 7 can be established 
by two methods (cf. A. T., p. 52, footnote) : 

(a) It is the difference between the resultant of the pressures 
on a flat plate, which is normal to the plate, and the total net 
force acting on a “reference surface” surrounding the plate, which 
is normal to the flow. This method is inconsistent (cf. A. T., 
p. 51, last paragraph) with the first order approximation used to 
determine the pressure distribution for a compressiole flow. In 
addition, a direct force determination on a reference surface, 
valid to the second order, does not appear to be available, so that 
this method will not be used hereafter. 

(b) The suction force is the limit of the net force acting on a 
small contour passing in the vicinity of the leading edge when this 
contour tends to coincide with the airfoil. From the first order 
approximation theory, pressure changes near a leading edge 
moving in an incompressible flow, or having a subsonic normal 
velocity component in a compressible flow, vary inversely as the 
square root of the distance to the edge. The pressure changes p 
have opposite values on the two faces of the plate, and there is no 
suction force, because of the antisymmetry of the p-distribution 
given by the linearized Bernoulli formula, 


p+ pUu =0 (3) 


where u is given by Eq. (2). Thus, the value of Eq. (1) for T is 
not the consequence of Eqs. (2) and (3), and a second order ap- 
proximation is required. [n this case, the correct Bernoulli 
formula, 

p+ + u)? — = 0 (4) 


replaces Eq. (3), and p contains a symmetric term —1*/2u%p, 
which shows that p is inversely proportional to the first power 
of the distance to the leading edge. In this case, there exists a 
finite suction force, as shown by R. Grammel in ‘Die hydro- 
dynamischen Grundlagen des Fluges’” (Braunschweig, 1917, p. 
21). The exact value of this force is Eq. (1). It is worth re- 


marking that Eqs. (3) and (4) give the same resultant of the pres- 
sures acting on both sides of a plate, which is equal to 2Uup. 


the Editor 


Thus, since Eq. (1) is the consequence of Eqs. (2) and (4), and 
not of Eqs. (2) and (3), the general expression given by the 
authors [Eq. (21)] must be checked for consistency with thesecond 
order approximations of all the quantities that appear init. This 
means that the variations of density with velocity must be taken 
into consideration. However, the actual equation of state must 
be replaced by an expression ‘compatible with the linearized dif- 
ferential equation of motion, which contains the term (1 — M%), 
where M is considered as a constant which appears in the equa- 
tion of state, which becomes 


dp/dp = U?/M? = Ct (5) 


In terms of the local velocity g, which has the components 
(U + u, v, w), Bernoulli equation is 


(dp/p) +qdq =0 (6) 


Eliminating dp between Eqs. (5) and (6) and integrating, it is 
found 


p/po = pa/ = exp. — (g? — U2) M?/2U? (7) 


where pq and po are absolute pressures. Thus, ~, is equal to 
zero when g is infinite, and the relative pressure has the finite 
value — pp at the leading edge. Thus, there is no suction force 
when a second order correction is introduced. 

Consequently, it is proved that, near the leading edge, second 
order corrections to the linear compressible fluid theory lead to 
finite pressures that are incompatible with the existence of suction 
forces. On the other hand, the first order theory, although lead- 
ing to infinite pressures along subsonic leading edges, does not 
imply the existence of suction forces, because of the antisymmetry 
of the pressure distribution. In addition, the strength of the 
singularity is insufficient to cause a finite force, since SI, *(C//z)dx 
tends toward zero with x. The complete theory cannot lead to 
infinite suction forces, even for perfect fluids, for the same reasons 
as the second order theory: the density decreases with pressure, 
and infinite negative pressures cannot be reached, contrary to 
the case of incompressible fluids. For this reason, d’Alembert 
paradox does not hold for compressible fluids, and the determina- 
tion of forces based on the application of this paradox may be 
incorrect. It is probable that the greatest errors occur in the 
case when there is a sharp edge, since in this case, the flow differs 
greatly from a potential flow. 

In conclusion, it appears established that, on the one hand, 
neither an approximate theory (first or second order) nor an ex- 
act one leads to the existence of suction forces along sharp sub- 
sonic leading edges; and, on the other hand, none of these theo- 
ries requires, for reason of consistency, the introduction of these 
forces. It remains to be seen, however, both experimentally and 
analytically, what relation exists between the sharpness of the 
leading edge and the drag of a delta wing. 


LEON BESKIN 
Design Research Division 
Bureau of Aeronautics 
Navy Department 
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Dependence of the Stress-Strain Curves of 
Cold-Worked Metals Upon the Testing 


Direction’ 


L. J. KLINGLER? anv G. SACHSt 
Case Institute of Technology 


ABSTRACT 


The yield strength and stress-strain curve of a cold-worked 
aluminum-alloy plate depend considerably upon the angle 
between the directions of the largest principal strains in cold- 
working and testing, respectively. This directionality develops 
as a result of a small amount of cold-work and then remains 
approximately constant. The effect is closely related to the 
Bauschinger effect. 


INTRODUCTION 


i bes RELATIONS BETWEEN STRESSES AND STRAINS for 
an annealed and presumably isotropic metal can 
be expressed as a universal function that is independent 
both of the direction of the test specimen and the state 
of stress (within certain limits). On the contrary, the 
stress-strain curve of metal, cold-worked by even a 
small amount, has been found to vary considerably 
in subsequent testing, depending upon the differences 
between either the directions or the states of stress in 
cold-working and: testing, respectively. 

The variations in the stress-strain curves are ap- 
parent primarily in the range of small plastic strains, 
up to approximately 2 per cent, as illustrated by 
Figs. 1 and 2 for 24ST and 24SRT aluminum-alloy 
sheet.! 

The dependence of the stress-strain curve of a cold- 
worked metal upon the direction of testing has ap- 
parently not yet been systematically investigated. 
An indication of this effect is contained in an ANC 
Bulletin.? According to Table I-1 of this bulletin, 
sheets in the commercial aluminum alloys 24ST and 
24SRT possess practically equal tensile strengths in 
the longitudinal and transverse directions. However, 
their yield strength (0.2 per cent) in tension is specified 
approximately 15 per cent lower in the transverse than 
in the longitudinal (rolling) direction, this fact being 
obviously based on the average of a large number of 
tests. The above-mentioned table shows, further- 
more, that the yield strength in compression of the 
sheet metals is approximately 12 per cent larger in 
the transverse than in the longitudinal direction. 


Received September 24, 1947. 

* The work reported in this paper was made possible by the 
Frederick Cottrell grant received from the Research Corporation, 
New York. 

¢ Research Laboratory for Mechanical Metallurgy, Case 
Institute of Technology (formerly Case School of Applied 
Science). 
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Fic. 1. Typical stress-strain curves for 24ST sheet. 
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Fic. 2. Typical stress-strain curves for 24SRT sheet. 
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Fic. 3. Stress-strain curves for cold-worked brass tested in 
tension. 


DISCUSSION 


The typical stress-strain curves in tension and 
compression for 24ST and 24SRT sheet, in Figs. 1 and 
2, clearly illustrate these differences in yield strength. 
The source of this anisotropy is the cold-working to 
which 24ST and 24SRT sheet are subjected as the last 
commercial operation. To obtain flatness, 24ST sheet 
is stretched and/or rolled by approximately 1 per cent, 
while 24SRT is rolled by 5 to 6 per cent to increase 
its strength. If no cold-work is applied, the stress- 
strain curves in both directions have been found to 
be almost identical in shape.’ In addition, aging at 
temperatures as low as 300° to 350°F. reduces the 
yield-strength difference considerably.' 

The temperature dependence of this anisotropy 
suggests that it is correlated with the development of 
residual stresses of the ‘second order.’’ * 

Such residual stresses are generally considered re- 
sponsible for the much-discussed differences in the 
stress-strain curves and, consequently, the differences 
in yield strengths of a metal prestrained in tension 
and tested in tension and compression (different stress 
states) parallel to the prestrain load. For a metal 
prestrained in compression and again tested in tension 
and compression parallel to the prestrain load, the 
effect generally known as “Bauschinger effect’’ is 
reversed. Figs. 1 and 2 illustrate this effect for 24ST 
and 24SRT, and Fig. 3 for brass. This latter graph 

* Residual stresses of the second order (microstresses) cannot 


be demonstrated mechanically as can residual stresses of the first 
order. They can, however, be shown to exist by means of 


X-rays.?-4 
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also shows that the stress relieving at a suitable tem. 
perature may eliminate the Bauschinger effect without 
affecting the hardness, the resulting stress-strain curve 
then being intermediate between those for the speci- 
mens prestrained by tension and compression, rte. 
spectively. 

Figs. 1 and 2 indicate, furthermore, that, after pre- 
straining by rolling, the transverse stress-strain curve 
in tension is similar to the longitudinal stress-strain 
curve in compression. Apparently, rolling acts, in 
this respect, similar to stretching by tension. 

The peculiar similarities between the anisotropy of 
the stress-strain curve in a particular stress state, and 
the dependence of the stress-strain curve upon the 
sign of the stress (Bauschinger effect) become more 
lucid if the strain state rather than the stress state is 
considered. The two phenomena can then be sum- 
marized by the statement that the stresses required for 
small plastic strains are lower the more the strain 
state in testing differs from that of the preceding cold- 
work. Previotls experimental evidence indicates that 
the Bauschinger effect first increases with strain and 
then remains approximately constant for strains 
exceeding 1 to 2 per cent, as illustrated in Fig. 4 for 
some tests, from various sources, on 24ST sheet pre- 
stretched in tension.! Fig. 5 shows, furthermore, 
that the anisotropy effect (i.e., the difference in the 
tensile yield strengths parallel and perpendicular to the 
direction of prestretching) depends upon the strain in 
much the same manner as the Bauschinger effect. 
To evaluate quantitatively the two different effects, 
it must be considered that the commercial sheet that 
was used was slightly anisotropic and also possessed 
slightly different values of yield strength in tension and 
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Fic. 4. Effect of prestrain in tension on compressive and tensile 
yield strengths in longitudinal direction for 24ST clad. 
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STRESS-STRAIN CURVE OF COLD-WORKED METALS 


compression, respectively. Therefore, if the four 
curves in Fig. 5 are shifted to possess the same initial 
yield strength, it can be seen (Fig. 6) that the magni- 
tude of the two effects discussed is practically identical 
for any given amount of*prestretching. 

In this paper, a few tests are reported which establish 
the anisotropy of the stress-strain relations for a metal 
prestrained in uniaxial tension. The aluminum alloy 
24SO, stretched by approximately 12 per cent in 
tension, has been selected in order to illustrate that the 
anisotropy is not restricted to 24ST but that it is a 
universal phenomenon. According to Figs. 4 to 6, 
the magnitude of the effect should depend little upon 
the amount of cold-work, provided that 1 per cent 
cold-work is exceeded.* 


RESULTS 


The various experimentally obtained stress-strain 
curves are represented by the points in Fig. 7. At 
small strains the stress values vary considerably with 
the testing direction, as was expected. For any strain 
larger than approximately 1.5 per cent, however, the 
observed stress values were independent of the testing 
direction within the limits of scattering, yielding a 
stress of 36,000 = 1,000 Ibs. per sq.in. at a total strain 
of 2 per cent. 

To correct for the small but undesirable experimental 
variations, the trend curves in Fig. 7 were not drawn 
through the test points but to conform exactly to a 


* The effects of large strains that result in noticeable crystallo- 
graphic directionality are not considered here. Such effects are 
superimposed to those discussed above. 
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Fic. 5. Effect of prestrain in tension on compressive and 
tensile yield strengths in longitudinal and transverse directions 
for 24ST (heat-treated from 24SO clad). 
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Fic. 6. Effect of prestrain in tension on compressive and 
tensile yield strengths curves from Fig. 5 corrected for initial 


anisotropy. 
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Fic. 7. Stress-strain curves in tension at various’ angles to the 
prestraining direction for 24SO prestrainedfin tension. 
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Fic. 8. Stress-plastic-strain curves in tension at various angles 
to the prestraining direction for 24SO prestrained in tension. 


stress of 36,200 Ibs. per sq.in. at 2 per cent strain. 
Eventually, the elastic strains have been deducted from 
these curves to yield stress-plastic-strain curves for 
each direction tested, 

These (corrected) stress-strain curves are assembled 
in Fig. 8. In addition, in Fig. 9, the stresses for various 
plastic strains, 0.05, 0.1, 0.2, 0.5, 1.0, and 2.0 per cent, 
are plotted as functions of the direction angle—i.e., 
the angle between the directions of stretching and 
testing. For any given strain the stress decreases 
first approximately linearly with the direction angle, 
to about 50°, and then decreases at a much slower 
rate, to 90°. The total difference increases rapidly 
with decreasing strain, approaching a value of 6,500 
Ibs. per sq.in. at zero plastic strain, or 24 per cent of 
the yield strength of longitudinal specimens. 

If a specimen is prestretched in tension by 12.0 
per cent, its transverse dimensions—i.e., at a direction 
angle of 90°—are reduced by 5.5 per cent. Conse- 
quently, any fiber located at direction angles over 
52.2° have been compressed, while those at smaller 
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Fic. 9. Effect of the angle between prestraining and testing on 
the stress required for various plastic strains. 


direction angles have been stretched. The results 
represented in Fig. 9 indicate that for a given plastic 
strain the stress decreases considerably with increasing 
direction angle until the respective prestrain becomes 
zero, and then changes only slightly. 

This statement is also in agreement with the fact, 
not discussed so far, that, after prestretching, the 
anisotropy in compression is considerably smaller than 
that in tension (see Figs. 1, 2, and 6). 
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A NenitaslGem hical Method of 
Characteristics for Anially Symmetric 
Isentropic Flow’ 


LESTER L. CRONVICHt 
The Johns Hopkins Unwersity, Applied Physi¢s Laboratory 


SUMMARY 


A graphical-numerical method of characteristics, based on the 
Tollmien! method, is developed for axially symmetric isentropic 
flow problems. Complete rules of procedure are presented for 
the application of the method. Similar methods have been given 
by, Ferrari? and Ferri.’ As an illustration, the flow around a 
conical-nosed body is calculated. In this problem a procedure is 
given for approximating the curved shock wave, assuming con- 
stant entropy behind the shock wave. 


INTRODUCTION 


METHOD for calculating steady, axially 
symmetric, supersonic flows was developed in 
1934 by Frankl.4 The method was an application of 
the Prandtl-Busemann® graphical method to problems 
with axial symmetry. Both the physical plane and the 
hodograph plane are used in the method. More re- 
cently, Sauer® has presented a numerical-graphical 
method that also uses both planes. 

At thé Fifth Volta Congress in Rome, 1935, Ferrari 
briefly outlined a method of characteristics for the solu- 
tion of such problems in the physical (or flow) plane 
alone. This method was given in more detail in later 
articles.2, Tollmien and Schafer! outlined a ‘‘field” 
method of characteristics in 1940. This paper outlines 
a “lattice-point’”’ method of characteristics which is 
essentially equivalent to that of Tollmien and Schafer. 
A similar method has also been given by Ferri.* The 
fundamental mathematical theory underlying the 
method of characteristics has already been sketched by 
the author in a previous article* and is not repeated 
here. 


DEVELOPMENT OF THE CHARACTERISTIC EQUATIONS 


For axially symmetric potential flow, the velocity 
potential ¢ satisfies the following differential equation 
in cylindrical r,z-coordinates :’ 


[1 — (u?/a?) — (2uv/a*)ore + [1 — (0?/a*) + 
(¢,/r) (1) 


where 


Received June 24, 1947. Revised and resubmitted November, 
1947. 

* This work was done at the Applied Physics Laboratory of 
The Johns Hopkins University, Silver Spring, Md., under Con- 
tract NOrd 7386 with the Bureau of Ordnance. 

Mathematician. 


= ¢, = z-component of velocity 

= ¢, = r-component of velocity 

= V[(e — 1)/2] — = 
sound 

q = Vu? + v*? = velocity 

Qm maximum velocity, corresponding toa = 0 

k ratio of specific heats (1.405 for air) 


ee 


local velocity of 


For supersonic flow, the two families of characteristics 
are solutions of 


[1 — (u?/a?) |dr? + (2uv/a*)dr dz + 
{1 — (v?/a*)|dz? = 0 (2) 


or 


dr/dz = (uv + aVq — a? — a*)/(u? — a’) (2a) 


Let @ be the inclination of the velocity vector to the 
axis of symmetry so that 


= 
(3) 
v = qsin@ 


Also let a be the Mach angle, defined by 
sin a = a/q (4) 


Substituting Eqs. (3) and (4) in Eq. (2a) and simpli- 
fying, it is found that 


dr/dz = (sin 8 cos 8 = sin a cos a)/(cos? @ — sin® a) 


(2b) 
After several transformations this can be reduced to 


dr/dz = tan (0 = a) (2c) 


Thus, the tangents to the characteristic curves of 
the first family form the angle a with the streamlines, 
and the tangents to the characteristic curves of the 
second family form the angle — a with the streamlines. 
Consequently, we identify them with the well-known 
“Jeft-running”’ and “right-running’’ Mach lines. 

Along these characteristics, Eq. (1), expressed in 
terms of g and 0, has a form readily adaptable to com- 


putation. By differentiation, 
du = + (5) 
dv = + ddr | 

from which 
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1 


Fic. 1. Location of unknown point 3 from known points 1 and 2. 


$222 = du — (6) 
= dv — ¢$,,dz 

Multiply Eq. (1) by dr dz and replace ¢,,dz and ¢,dr 

by the values given in Eqs. (6). Collecting terms, this 

equation becomes 


(1 — u?/a*)du dr + (1 — v?/a?)dv dz + 
(v/r)dr dz — o2[(1 — u?/a?)dr? + (2uv/a*)dr dz + 
(1 — v?/a?) dz?] = 0 (la) 
But along the characteristics, the last term of Eq. 


(la) vanishes by virtue of Eq. (2). Thus Eq. (1a) be- 
comes 


(1 — u?/a*)(dr/dz)du + (1 — v?/a?)dv + (v/r)dr = 0 


(1b) 
Differentiating Eqs. (3), 
du = dq cos 6 — sin (7) 
dv = dq sin 60 + qd0 cos 0 


Substituting Eqs. (3), (4), (2c), and (7) in Eq. 
(1b) and simplifying, the following equations are 
obtained : 


cot a sin 6 sin a dr 


sin (@ + a) r (le) 
along 
dr/dz = tan (@ + a) 
and 
cot a sin 6 sin a dr 
along 
dr/dz = tan (0 — a) 
Let 
q q a 
where q’ = Also let 
F = sin @ sin a/sin (6 + a) 
G = sin @ sin a/sin (6 — oH (9) 
Then Eggs. (1c) and (1d) can be written 
dé — Qdq’ + (F/r)dr = 0 (le) 
dé + Qdq’ — (G/r)dr = 0 (1f) 


along the “‘left-running”’ and ‘“‘right-running”’ character- 
istics, respectively. 

It is observed that Q may be calculated beforehand 
as a function of g’. As an aid in the later computa- 
tions, curves of Q, a, and M may be plotted as functions 
of q’ from the following formulas, which can be readily 
derived from the equation defining a, below Eq. (1). 


M = V [2/(k — — 9’)] 
a are sin (1/M) (4a) 
Q = VM (8a) 


A sample of these curves is given in Figs. 2a and 2b. 


LATTICE-POINT METHOD OF SOLUTION 


The solution of Eqs. (1c) and (1d) is approximated by replacing the differentials by differences. Let points 1 
and 2 in the adjoining figure be known points—that is, their locations and values of g’ and @ are known. The new 
point 3 is located at the intersection of the ‘‘left-running”’ characteristic from point 1 and the ‘‘right-running”’ 
characteristic from point 2. The characteristics are approximated by line segments as indicated. Then, from 


Eqs. (le) and (1f), 


03 — — Qilgs’ — + (Fi/n)(rs — = (11) 
— + Qe(gs’ — ge’) — (Go/r2)(rs — 72) = OS. 


Solving Eqs. 11) for gs’, 


qs’ = + +h 


Either of the following formulas may then be used for 63: 


rs 
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(1¢e) 


(1d) 


(8) 


AXIALLY SYMMETRIC ISENTROPIC FLOW 


= 0, + Qilgs’ — — Fillrs — 


Preferably one formula is used as a check on the other. 


If point 1 lies on the axis of symmetry, then the last 
term of the first of Eqs. (11) becomes indeterminate, 
since 6; = 0 and 7, = 0. However, it may be readily 


evaluated. 


sin 6; sin Q — 


F 
fim — — = lim 


n—>o sin (0; r; 
o 
sind 63 
= £3 lim = 73 \ — = £3 — = 63 
dr/; r3 


In this case Eq. (13) becomes 
9 = Qs(gs’ — qu’)/2 (15) 
and Eq. (12) becomes 


+ Qe + 2G2 + 26, 
(16) 
+ 20, 


If a point lies on a fixed boundary, 63 is known and 4;’ 
can be determined from the appropriate equation of 
Eqs. (13) and (14). In this connection, the axis of 
symmetry can be considered as a fixed boundary for 
which 6 = 0. 

In the case of jet flow into air at lower pressure, the 
external pressure ratio p/p» is given, where fp is stag- 
nation pressure. Thus, the value of Mach Number 
and consequently of g’ is known on the jet boundary. 
The boundary direction can then be obtained from the 
appropriate equation of Eqs. (13) and (14). 


RULES OF PROCEDURE 


The basic rules of procedure are given for obtaining 
anew point 3 at the intersection of characteristics from 
two known points 1 and 2. It is assumed that the 
flow variables g’ and @ are known at the given points. 
All computations and construction are made for the 
flow pattern above the axis of symmetry. 

(1) From Figs. 2a or 2b determine a and ay corre- 
sponding to the given values g;’ and q,’, respectively. 
Draw a straight line segment from point 1, making an 
angle 6; + a, with the axis of symmetry, and another 
straight line segment from point 2, making an angle 6. — 
a With the axis of symmetry (Fig. 1). The intersection 
of these line segments locates point 3. 

Special case: Either point 1 or 2 may be near a 
boundary curve or axis of symmetry on which the new 
point 3 will be. In such a case only one line segment 
isneeded. (See illustration.) 


— Qe(gs’ — ge’) + Ge[(rs — r2)/re] 


(13) 
(14) 
Note that the angles above are in radians. 
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with velocity ratio q’. 


. Variation of Mach Number M, Mach angle a, and Q 
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Fic. 2b. Variation of Mach Number M, Mach angle a, and Q 


with velocity ratio q’. 
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Fic. 3. Pattern of characteristics for conical-nosed body. 


(2) Calculate gs’ angi 03 from Eqs. (12), (13), and (14) 
unless point 1 lies on the axis of symmetry, in which case 
Eqs. (14), (15), and (16) are used. 

(3) From Figs. 2a or 2b determine a;, Q;, and M, 
corresponding to q;’ found in Step 2. Also measure the 
distance r;. If the construction is carried out on rec- 
tangular coordinate paper, this distance may be read 
at once. 

(4) Iteration: If Mz; is low, it is usually necessary 
to follow a method of iteration in order to improve the 


accuracy in the determination of the flow pattern. Of © 


course, the application of the iterative procedure de- 
pends on the degree of accuracy desired. The following 
rules apply: 

(a) Calculate the average quantities: 


(ay + as)/2, (az + a3) /2 
(4, + 63)/2, 82, c= (0. + 63) /2 
(n+ r3)/2, 74, (ro + 13)/2 


(b) From Figs. 2a or 2b determine Q(a;,3) and 
Q(ae, 3), corresponding to a, 3 and ae, 3, respectively. 

(c) Calculate F(a, 3) and Glas, 3, %, 3). The 
quantities Q(a1, 3), Q(a2, 3), s)/1, 3, and G 
(a2, 3, 92, 3)/re, 3 represent values of the coefficients of 
dq’ and dr in Eqs. (le) and (1f) at intermediate points 
of the characteristic curves. These evaluations are 
usually more nearly correct than the values at either 
end point. 

(d) Draw new line segments from points 1 and 2 
making angles 6;, 3 + a1, 3 and 3 — as, s, respectively, 


with the axis of symmetry. The intersection of these 
lines locates a point 3’, which is more nearly coincident 
with the intersection of the curved characteristics from 
points land 2. Measure the distance r3’. 

(e) Using the applicable formulas of Step 2, calculate 
gs’ and from the known values 61, and and 
the intermediate coefficients found in (c) above. 

(f) Find ay, My, and Qs on Figs. 2a or 2b corre- 
sponding to q3’’. 

(g) If the location of point 3’ and the flow variables 
at 3’ do not differ appreciably from the location of 3 and 
its flow variables to the accuracy desired, the pro- 
cedure may be stopped. Otherwise the procedure is 
repeated to locate point 3”, etc. 

(5) It should be noted that the angles used in the 
preceding formulas are in radians. 


ILLUSTRATION 


Flow Over a Conical-Nosed Body 


To illustrate the procedure outlined in the preceding 
section, the characteristics pattern about a conical- 
nosed body is constructed (see Fig. 3). The body is 
immersed in a steady flow at initial Mach Number 
2.047. The conical nose has a semivertex angle of 
30°. At the point (3, 1.73) the conical nose becomes 
tangent to a circular arc of radius 10 in the meridian see- 
tion of the body. This arc may be used as a transition 
between the conical nose and a cylindrical body. 
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From integrations performed by the Taylor-Maccoll 
method,* the conical flow field is completely deter- 
mined. The set of integrations referred to in reference 
9 provides sufficiént data to serve as the initial values 
necessary for illustrating the method of characteristics. 
In particular, it is found that the original shock-wave 
angle is 47.4°. Furthermore, values of wu’, 0’, and a’? 
(radial and tangential velocity components, and square 
of sound velocity, respectively, referred to maximum 
velocity, Gm) are tabulated for various steps in the 
integrations. From these values the data of Table 1 
may be obtained with the use of some well-known 
formulas. The angle ¢ is the angle between the axis 
of symmetry and a radial line from the vertex of the 
cone. The data for e = 47.4° are given for the down- 
stream side of the shock wave. 


TABLE 1 
Pt q’ 6° a? (0 + clea.” M 
1 30 0.500 30.0 51.2 tay 1.280 
76.2 
e - 40 0.511 22.0 49.2 1.320 
68.6 
4 45 0.523 18.8 47.2 1.363 
64.65 
47.4 0.530 17.2 46.1 1.388 


The first task is the construction of the character- 
istic at point 1 which separates the conical flow field 
from the rest of the flow field. Since the Mach Num- 
bers are low, the process of iteration is necessary in 
this problem in order to obtain a fair degree of accuracy. 
More properly, a finer mesh of characteristics should 
also be used. However, in order to show more clearly 
the graphical construction (including some of the itera- 
tion process), the mesh has purposely been kept large. 
Only the first and last approximations are shown in 
Fig. 3. As an approximation, the line segment drawn 
from point 1 to point 2 is drawn at an angle 


(@+ a),2= (6+ (6+ a)2]/2 


ie., at an average (0 + a) between the end points. The 
same procedure is followed for the other line segments 
in the characteristic 1-2-4-7._ The angles are given in 
Table 1. 

The flow field in the triangular region 1-7—10 is isen- 
tropic and is calculated by the method of iteration 
described previously. The iterative procedure is 
stopped at any point (in this problem) if the new value 
of g’ differs from the preceding value by no more than 
0.002 and if the new value of @ differs from the pre- 
ceding value by no more than 0.1°. This choice, of 
course, is always dependent on the accuracy desired. 
Values of 6 along the boundary can be calculated or 
measured graphically. The numerical calculations for 
the region 1-7-10 are carried out conveniently by 
means of a tabular form as shown in the first part of 
Table 2. 


APPROXIMATE CONSTRUCTION OF CURVED SHOCK WAVE 


At point 7 the shock wave curves because of the 
effect of thg change in slope of the boundary. Curva- 
ture of the shock wave introduces vorticity for which 
the method of this article is not valid. However, for 
free-stream Mach Numbers that are not too high and 
bodies that are not too sharply curved, a good approxi- 
mation to the true solution can be obtained by assum- 
ing that the flow continues irrotational downstream of 
the shock wave. This is equivalent to assuming uni- 
form entropy changes across the entire length of the 
shock wave, so that the flow behind the shock wave is 
isentropic and the method of characteristics given 
above can be used for the solution. 

The next task is the construction of the shock wave 
beyond its straight portion. Using straight-line inter- 
polation choose an intermediate point 7a on the “‘right- 
running’’ characteristic 7-8’. The location of point 
7a will determine the size of the line segment of shock 
wave drawn at point 7 to a new point 11. Basically, 
a point 11 must be located on the “‘left-running”’ char- 
acteristic from point 7a, and on the shock line from 
point 7, such that the flow variables calculated by 
shock wave theory (with given free-stream Mach 
Number) and characteristics theory (from known point 
7a) will agree to the accuracy desired. This can be 
accomplished by an iterative procedure in the following 
steps: 

(1) Extend the original shock line through point 7; 
then draw the ‘‘left-running” characteristic from point 
7a using the angle (9 + a);, until it intersects the shock 
line at point 11. 

(2) Knowing free-stream Mach Number (2.047) and 
shock-wave angle (47.4°), the downstream flow vari- 
ables at point 11 can be determined from shock-wave 
tables'’: |! or graphs such as Figs. 4and 5. In Table 2, 
these values are given for point 11 to distinguish them 
from characteristics data for point 11. Note that the 
data for point 11 can be obtained from point 7 since the 
shock-wave angle is the same. 

(3) Assuming 6;, = 67 (i.e., flow direction obtained 
from shock-wave theory is assumed for point 11 in char- 
acteristics theory), calculate gu’ from characteristics 
theory and obtain the corresponding M from Fig. 2a or 
2b. 

(4) Compare M), with My. If the values differ, as- 
sume 


My = (Mu + Mji)/2 


as a new value of Mach Number after the shock wave, 
and determine the corresponding flow direction @j- 
from Fig. 5. 

(5) From point 7 draw a line segment at the shock- 
wave angle found in Step (4). This represents the next 
approximation to the shock line. 

(6) Assume 0, = Oj. Draw the “left-running”’ 
characteristic from point 7a at an angle 6,,. + a,,, where 
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Bay, = (74 + 655") /2 


and 
Gay, = + 


The intersection with the shock line drawn in Step (5) 
jocates point 11’. 

(7) Calculate g,, from charactéristics theory, using 
coefficients evaluated at an intermediate point (indi- 
cated by 7a-11’ in Table 2), and then obtain My, from 
Fig. 2a or 2b. If this value does not agree with the 
assumed value of Mj, the procedure is repeated from 
Step (4). 

The flow field between the “right-running”’ charac- 
teristics from points 7 and 11’ can be calculated ac- 
cording to the method described for the region 1—7-10. 
An illustration is given for point 12. Other segments 
of the shock line such as 11’—16’ may be drawn accord- 
ing to the procedure outlined above. 


COMMENTS ON THE METHOD 


It is observed that the quantity sin (@ + a) appears 
in the denominator of the third term of Eq. (lc), and 
sin (9 — a) appears in a similar position in Eq. (1d). 
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Fic. 4. Mach Number M;, before shock vs. Mach Number M; 
after shock for various flow deflection angles 8. 


This might appear to lead to a singularity in each equa- 
tion under certain circumstances. For example, if 
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- 6 = —a in Eq. (1c), the denominator vanishes. But 


this requires that the ‘‘left-running’’ characteristic, 
which makes an angle (9 + a) with the axis of sym- 
metry, be parallel to the axis of symmetry. Hence, 
dr = 0 and, in reality, the term becomes indeterminate. 


In this exceptional case it may be shown that the term 


approaches a limit given by — [(sin® a)/r]dz. Similarly, 
if @ = a in Eq. (1d), the third term also approaches a 
limit given by — [(sin* a)/r]dz. Since these phenomena 
appear so infrequently, no special calculation formulas 
are listed for these cases. 

The method can be readily applied to the design of 
axially symmetric supersonic nozzles following a pro- 
cedure fundamentally similar to that of Frankl.‘ That 
is, the transition flow field between a three-dimensional 
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source flow and a uniform parallel flow may be calecy. 
lated by the ‘‘lattice-point’’ method of characteristics, 
Any streamline of this field may then be used as q 
boundary for directing an initially radial flow (the 
source flow) into a parallel flow. The drawing of the 
streamlines requires some method of interpolation be- 
tween the calculated points on the flow field. If a 
sufficient number of points are calculated, linear inter. 
polation will suffice. 

Although the procedure above has been given for 
irrotational flows, the method can be applied to flows 
with vorticity. This application has been given 
by Ferri.? It requires an assumption as to the 
distribution of vorticity downstream of the shock 
wave. 
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Thermal Fin Effect in Heat Anti-Icing 
Corrugations 


JULIUS JONAS* 
Northrop Airéraft, Inc. 


SUMMARY 


The importance of heat transfer detail design in heat anti-icing 
installations is stressed with particular reference to the design of 
conventional inner skin corrugations. An analysis is made in an 
effort to evaluate the thermal fin effect contribution made by the 
inner skin corrugations. The results obtained emphasize the 
need for careful detail design of conventional inner skin corruga- 
tions so as to obtain maximum anti-icing potentialities for a given 
installation. 


NOMENCLATURE 


a = distance AB in Fig. 1, ft. 

b = distance ABC in Fig. 1, ft. 

Ci, C2, C3, Cs, Cs, Ce = integration constants 
d = distance ABD in Fig. 1, ft. 


fi = corrugation internal heat transfer coefficient, B.t.u. per 
hr. per sq.ft. per °F. 

f. = effective external heat transfer coefficient based on the 
total heat dissipated, B.t.u. per hr. per sq.ft. 
per °F. 


fet. = effective internal heat transfer coefficient [defined by 
Eq. (34)], B.t.u. per hr. per sq.ft. per °F. 

k = skin conductivity, B.t.u. per hr. per sq.ft. per °F., ft. 

Ki = constant defined by Eq. (13), ft.~! 

K. = constant defined by Eq. (14), ft.~! 

K; = constant defined by Eq. (15), °F. 

K, = constant defined by Eq. (16), ft.~! 

m = constants defined by Eq. (26a) through Eq. (26i) 


Qz = heat conducted in skin along x-axis, B.t.u. per hr. per 
sq.ft. 

Qz, = heat conducted in skin along x-axis, B.t.u. per hr. per 
sq.ft. 

Q:, = heat conducted in skin along x-axis, B.t.u. per hr. per 
sq.ft. 

R = fes./f; = heat transfer effectiveness 

t; = corrugation heated air temperature, °F. 

t, = free-stream temperature, °F. 

tz = skin temperature along x-axis, °F. 

= skin temperature along x-axis, 

tz, = skin temperature along x-axis, °F. 

t, = average external skin temperature defined by Eq. (33), 

x = distance variable (see Fig. 1) with origin at point “A,” 
extending from point “A” to point “B” 

x; = distance variable (see Fig. 1) with origin at point ‘‘A,”’ 
extending along ABC 

x2 = distance variable (see Fig. 1) with origin at point ‘‘A,”’ 
extending along ABD 

6; = corrugation skin thickness, ft. 

62 = outer skin thickness, ft. 


6s = + ft. 


Subscript 
av. = average 


Received June 12, 1947. 
* Aeronautical Engineer. 


INTRODUCTION 


iy DESIGNING conventional inner skin corrugations 
for surface heat anti-icing installations, attention 
appears to have been given in most designs only to 
the selection of the size of the gap existing between the 
surface outer skin and the inner skin corrugation (see 
Fig. 1). The effects of width, pitch, and thickness of 
the inner skin corrugation on the heat transfer have 
hitherto been neglected. 

With heat anti-icing requirements assuming ever 
greater importance, it becomes essential to improve on 
heat transfer detail design so as to attain maximum 
anti-icing potentialities. The selection of the corruga- 
tion gap size will primarily establish the internal heat 
transfer coefficient pertaining to the heat flow from the 
heated air to the surfaces enclosing the passage. 

“The heat flow may be pictured as follows: Heat is 
transferred rom the heated air to the outer skin di- 
rectly. It is further transferred from the heated air to 
the inner skin corrugation, whence it is conducted into 
the outer skin. In that respect the inner skin corruga- 
tion provides a fin effect, effectively increasing the heat 
transferred from the heated air to the outer skin. This 
will result in increased outer skin temperatures, which 
will permit the handling of icing conditions of increased 
severity. Thus, in designing the heat anti-icing corru- 
gation passages, it appears that particular attention 
should be given to design details affecting the magni- 
tude of the inner skin fin effect. 

A heat transfer analysis is made below, correlating 
the parameters affecting the magnitude of the thermal 
fin effect of the inner skin corrugations. A ‘“‘heat 


<—— ASSUMED PERFECTLY INSULATED 
INNER SKIN CORRUGATION——, 
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transfer effectiveness’’ is defined and used as a measure 
of the thermal fin effect. 


HEAT TRANSFER ANALYSIS 


In the following heat transfer analysis it is assumed 
* that no heat is conducted through the inner and outer 
skin chordwise. This assumption leads to a two- 
dimensional heat flow analysis in a spanwise plane (see 
Fig. 1). It is further assumed that the inner skin 
corrugation is perfectly insulated, as shown in Fig. 1. 
With some complication of the expressions involved, 
the analysis could readily be extended to account for 
transfer of heat across the surface, assumed here as 
perfectly insulated. There is, however, a practical 
justification for this assumption. 

This surface is primarily exposed to heated air flowing 
spanwise in a conventional supply duct formed by a 
shear web and the chordwise inner skin corrugations. 
The spanwise air-temperature drop in the supply duct 
must be kept to a minimum so as to provide adequate 
anti-icing for the outboard surface regions. This in 
turn requires a minimum of heat transfer from the sup- 
ply duct air into the inner skin corrugations. 

It is further assumed that a perfect thermal bonding 
exists between the inner and outer skin. 

The heat conducted through a skin increment dx is 
given by 


(Qz)o (k/dx)6,dt, (1) 
(Qz,)a cine = > (2) 
(Ox)a cre ca = —(k/dx,)b3dtx, (3) 


Differentiating Eqs. (1), (2), and (3) with respect to x, 


(dQ,/dx)o — kb; (d*t,/dx?) (4) 
(dQx,/dx1)q = — (5) 
(dQx,/dx2)q <a = (6) 

(dQ,/dx) values may also be written: 
(dQ,/dx)o —filte (7) 


(dQx,/dx1)q —filtz, ty) felte, te) (8) 
(dQz,/dx2)q —feltxy te) (9) 


Equating Eqs. (4) and (7), Eqs. (5) and (8), and Eggs. 
(6) and (9), respectively, and rearranging, 


(d*t,/dx*) — (f;/R6,)te + = 0 
d*tr, fi +S, Sit: 
dx,? ( hbo Vie + kb» = 0 (11) 
— (f./Rds)tr + (fe/kbs)te = 0 (12) 
Indicating by 
= Ky? (13) 
(fe + fe)/Rb2 = (14) 
(fits + Sete) /(fi + fe) Ka (15) 
f./kis = Ke? (16) 


Differential Eqs. (10), (11), and (12) are readily solved: 
<e™ Ce" + Coe + ty 


(17) 
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Il 


= + Ce + K, 


The six coefficients in Eqs. (17), (18), and (19) may be 
reduced to three by the following three boundary condj- 
tions: 


(18) 
(19) 


(dt,/ dx), = 0 (20) 
(dtz,/dx;)2,-» = 0 (21) 
(dtz,/dx2)z,-4 = 0 (22) 


Solving Eq. (17) for boundary condition [Eq. (20)], Eq. 
(18) for boundary condition [Eq. (21)], and Eq. (19) 
for boundary condition [Eq. (22)], it readily follows 
that 


= C 
CG, = 
Ce 


Three additional boundary conditions are required to 
solve for Ci, C;3,and C;. The following boundary condi- 
tions are physically justified at x = x; = x2. = a: 


= (23) 


Introducing the following constants, 


my, = + (26a) 
mrs = 4. (26b) 
myo = K; — t, (26c) 
My = My (26d) 
Mes = eX — a) (26e) 
= ty — (26f) 
m3; = = (26g) 
M33 = eXs(2b—2)) (26h) 
ms = — (26i) 


C;, C3, and C; are readily determined by solving Eqs. 
(23), (24), and (25) simultaneously: 


+ M291113M 35 


= (27) 
+ — 1311113725 
C. MyMo9M35 — + M3130 (28) 
3 => 
+ M2 — 11131111325 
— + + 113111131129 
C; = (29) 


Having thus determined all constants in Eqs. (17), (18), 
and (19), respectively, it is now possible to compute 
average skin temperature values as follows: 


(trav.)o “ete aK (e + ty (30) 
1 
Cs K:(2b—a) 
<4 <0 = (b — (e (31) 


Cs 


(32) 
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THERMAL FIN EFFECT IN HEAT ANTI-ICING CORRUGATIONS 


Finally, an average external skin temperature value /, 
over portion DBC (see Fig. 1) is given by 

(b — a) Ja <n<ot(d— a) (tr, a 
b+d—2a 
Having thus established the temperature pattern of the 


(33) 


Analytically, ’ 
fea (b + d — 2a)(t, — t.) = falt, — 


and finally introducing the ratio 


o<e<al + fib — — 


alt, x (tray. o<z<al + (b ad a)(t, 


165 


inner and outer skin, it is possible to define an effective 
internal heat transfer coefficient f.g as follows: fem is 
taken as a heat transfer coefficient that would be related 
to the heat transferred by the hot air of temperature ¢; 
to an outer skin portion DBC = 6 + d — 2a in order to 
maintain it at a temperature /, as defined by Eq. (33). 


(tray. a << o] 


(34) 


This ratio R is termed a “heat transfer effectiveness”’ 
ratio and is a measure of the thermal fin effect contribu- 
tion made by the inner skin corrugation. 


RESULTS 


In order to discuss the results thus obtained, several 
particular cases were investigated. For all cases a con- 
stant gap of 0.125 in. was used. 


(1) Effect of variation of (6 — a) and (d — an on heat 
transfer effectiveness R. (See Fig. 2.) 


As shown in Fig. 1, (5 — a) = BC represents the half 
width of the corrugation passage, and (¢d — a) = BD 
represents the half width of the space between two 
consecutive corrugations. Fig. 2 permits the following 
conclusions: For BD = 0 the heat transfer effective- 
ness increases with a decrease of BC and approaches 
infinity at BC = 0. But this is clearly a trivial case. 
In most installations BD # 0 and is big enough to 
provide a rivet attachment to the outer skin. BC values 
have been ranging anywhere from 0.375 to 1.00 in. and 
more. The larger values were often selected in an effort 


to keep the production cost down. 


Fig. 2 clearly indicates that for a given value of BD 
there is a small range of BC values for which the heat 
transfer effectiveness becomes a maximum. The ex- 
istence of a maximum value of the heat transfer effec- 
tiveness at a given value of BC may be explained as 
follows: For small values of BC a large portion of the 
heat conducted from the inner skin must be used for 
the heating of the portion BD. With increasing values 
of BC larger and larger proportions of the heat con- 
ducted from the inner skin into the outer skin can be 
used effectively to raise the temperature of the portion 
BC of the external skin. This would indicate that the 


heat transfer effectiveness would continue increasing 
with an increase of BC. Actually, however, with in- 
creasing values of BC less and less heat conducted into 
the outer skin manages to reach the regions near point 
“C,”’ because it is being lost to the outside before it can 


(6 + d — 2a)(t, — 
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thermal fin effectiveness. 


THERMAL FIN EFFECTIVENESS — R 


reach point “C.” Thus, at a certain value of BC, a 
maximum value of heat transfer effectiveness is reached 
and beyond this a decrease is experienced. 

The use of BC values of less than the optimum value 
is obviously undesirable. By the same token, values of 


BC larger than the optimum value are undesirable not 
only because it results in a loss of heat transfer effective. 
ness but also because it results in the building up of 
higher and higher temperatures in the region of point 
‘“‘A”’ (see Fig. 1) which increases the buckling potential 
of the inner skin. Thus it appears desirable to keep 
BD values to a minimum and to select the correspond. 
ing optimum BC value. 

(2) Effect of inner and outer skin thickness on the 
heat transfer effectiveness. (See Fig. 3.) 

Fig. 3 indicates that the inner skin thickness is more 
critical than the outer skin thickness. However, over 
the range of practical inner and outer skin thickness 
values, no appreciable improvement in heat transfer 
effectiveness is obtained through an increase of skin 
thickness which may be warranted by a corresponding 
weight increase. 

(3) Effect of internal and effective external heat 
transfer coefficient on the heat transfer effectiveness, 
(See Fig. 4.) 

Fig. 4 indicates that over a useful range of heat trans. 
fer coefficients the heat transfer effectiveness will 
slightly decrease with a decrease of f,/f; for a given 
value of f,. It will decrease with an increase of f, fora 
given value of f./f;. Generally speaking, as the heat 
required to'maintain a given skin temperature increases, 
the heat transfer effectiveness will decrease. 


Letter to the Editor 


Dear Sir: 


In the Letter to the Editor by Professor Krzywoblocki and me, 
which appeared in the October, 1947, issue of the JOURNAL, two 
errors have since been determined. On page 564 in the second 
line below Eq. (1), the definition of V, should have read “‘ve- 
locity of propeller slipstream, with respect to the propelled body, 
ft. per sec.” . 


The second error appears on page 566 in Table 1. The thrust 
at maximum mp, for a rocket should read “‘(W;/g)v;”’ instead of 
“(W7/g) V.” 


R. W. McCroy 
Associate Professor 
Department of Aeronautical Engineering 
University of Illinois 
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SUMMARY | 


The flow at the rear of a two-dimensional airfoil moving at 
supersonic speed is of interest from a stability standpoint, 
since the tail surfaces are likely to be influenced by this flow 


field. 
The flow angle immediately behind the trailing edge of a flat 


plate airfoil and an airfoil with thickness has been calculated 
herein on the basis of oblique shock and Prandtl-Meyer wave 
theory. The deflection angle of the flat plate airfoil at a given 
Mach Number is shown to first order to vary as the fourth power 
of the angle of attack. The magnitude of the flow deflection 
between Mach Numbers of 1.2 to 5.0 is extremely small, being 
aways less than 0.06° for an angle of attack of 0.1 rad. 
(5.73°). 

'YThe mechanism of the return of the flow downstream of the air- 
foil to the free-stream direction is also discussed. 

The problem of the flow deflection immediately rearward of the 
two-dimensional supersonic airfoil has been previously treated 
in R. & M. No. 1930; the results of that paper, however, appear 
to be in error. 


INTRODUCTION 


ee FLOW AT THE REAR of a two-dimensional airfoil 
moving at supersonic speed is of interest from a 
stability standpoint, since the tail surfaces are likely 
to be influenced by this flow field. In calculations of 
the flow past two-dimensional supersonic airfoils, it is 
generally assumed that the flow just behind the trailing 
edge is in the free-stream direction. This assumption 
does not affect the surface pressures calculated but, as 
far as the wake direction is concerned, is equivalent 
to regarding the flow about the airfoil as completely 
isentropic. 

The flow angle immediately behind the trailing edge 
of a flat plate airfoil and an airfoil with thickness has 
been calculated herein on the basis of oblique shock and 
Prandtl-Meyer wave theory. The mechanism of the 
return of the flow downstream of the airfoil to the free- 
stream direction is also discussed. Although the effect 
of viscosity upon the flow angles is probably large, it is 
important to know the magnitude of the deflection for 
the ideal fluid. 


Received July 24, 1947. 
* Reproduction in whole or in part permitted for any purpose 
of the United States Government. This paper is an outgrowth 
of work done in supersonics at Princeton University for the Office 
of Naval Research, Project Squid. 
+ Research Associate, Department of Aeronautical Engineer- 
ing. 

t Assistant Professor, Department of Aeronautical Engineer- 
ing. 


The Flow at the Rear of a Two-Dimensional 
Supersonic Airfoil’ 


A. KAHANE} anp LESTER LEESt 


Princeton University 


The problem of the flow deflection immediately be- 
hind the trailing edge has been previously treated;! 
however, since the pressure changes through the oblique 
shock waves were considered as isentropic changes, the 
results of that paper appear to be in error. 


SYMBOLS 


= velocity of sound 

= stagnation velocity of sound 

= Mach Number, w/a 

static pressure 

= dynamic pressure, pw?/2 

= flow velocity 

= angle of attack . 

= flow deflection angle immediately rearward of trailing 
edge 


WR ese ee 


y = ratio of specific heats 

6 = wedge angle at leading and trailing edge of thick air- 
foil 

@ = flow direction angle 

p = mass density of fluid 


THE FLow DIRECTION IMMEDIATELY TO THE REAR OF 
THE TRAILING EDGE 


The Flat Plate Airfoil 


The flow immediately rearward of the trailing edge 
of a two-dimensional flat plate airfoil is first considered. 
If the airfoil is inclined to the freé stream at the angle 
of attack a, the change of flow direction at the leading 
and trailing edges is accomplished by two centered 
Prandtl-Meyer expansion waves and two oblique shock 
waves (see Fig. 1). The condition that must be satis- 
fied immediately rearward of the trailing edge is that 
the fluid in regions 4 and 5, which are in general sepa- 
rated by a vortex line AA’ of discontinuity of velocity, 


Fic. 1. The flat plate airfoil. 
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“a density, temperature, and entropy, must have the same Ap = a(+C0 + Ce? + C,0°'+ C64 = ...) 
flow direction and static pressure. 
The pressure change through isentropic waves was_ where q is the dynamic pressure ahead of the wave, 
expressed as a power series in the angle of deflection by ,w,"/2, and @ is the angle of deflection, positive in q 
Busemann? as counterclockwise sense. 


The signs to be used depend on the specific orientation of the surface and wave considered. For our purpose, 
a see the upper sign is used for the upper surface and the lower sign for the lower surface. The coefficients* C, are 
ee functions of the Mach Number ahead of the wave and are as follows: 


C, = 2/V mM? 
C, = [(M? — 2)? + yM*]/2(M? — 1)? 


ul? 4 _5+77- + + 1ly? — 8y — 
6 2(7 + 1) 24(y + 1) 4 — */s)? 
(Mm? — 1)” (M2 — 1)’ 


and the value of C;, is not needed here. 
Similarly, the pressure change through an oblique shock was expressed by Busemann as 
Ap = + C6? = (C, — + (C, — =...) 
where Dt is 


and E will not be needed. Thus it can be seen that the pressure change for isentropic compression and oblique 
shock waves differs only by terms of 6* and higher. 

Similarly, the change in dynamic pressure through isentropic and oblique shock waves can be expressed as a 
power series in deflection angle as: isentropic wave 


Ag = + = G0? + Got + ...) (1) 
oblique shock wave 
Aq = + G6? = (G, — + (G, — + ...] (2) 
where 
= (M? — 2)/V 1 
G. = {[(3 — y)/4]M* — 3M* + 4M? — 2}/(M? — 1)? 
6 —3y? + 177 — 36 —4y + 21 
an — 1)" — 1)” 
and G, and K will not be needed. 


* It should be noted that the values of C, and D as given by Busemann?® and Sauer? are incorrect. 
t See previous footnote. 


The flow direction behind the trailing edge, 8, will and 
now be expressed to terms of order as high as a‘. If all ps* — pi = u(—CiB + C6? — C38? + CyB4 — ...) 
the waves were isentropic, then 


where an asterisk over a symbol indicates its value for 
— pi = + + + C84 +...) a completely isentropic flow. 
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TWO-DIMENSIONAL SUPERSONIC AIRFOIL 


Since shock waves are not isentropic, it is evident 
that* 
= —@[D’(a + + E(a + +...) 


and 
ps — = + Eat+ ...) 


If the leading edge oblique shock on the lower surface 
were isentropic, the pressure difference in the expansion 
from 4 to 5 would be 


— = + B) + G"*(a + 8)? — 
+ + Ci"*(a + ...] 


whereas, actually, the pressure difference is 


ts— Ps = + B) + C"(a + — 
— + + (C” — + — ...] 


Expressing the fact that p; = ps or 


(bs — Ps*) + — pi) = (bs — fa) + (da — 
— pu*) + — pi) 


and substituting in the above equation gives 


2g + 68° + ...J = (a — 
+ mE'(a + B)* — + 
— quCi")(a + B) + — X 
(a + 8)? + — —D")}(a + + 
[qs*C."* — — (a + +... (3) 


It can be shown utilizing Eqs. (1) and (2) and by 
expanding C,” and C,”* in a Taylor’s series that 


— + B) = — + + 
O(a) + ...] = + B)a® + 0(a°) 


and that the terms to the right of that expression in Eq. 
(3) are of order a® and higher. Similarly, if D’ and E’ 
are expanded in a Taylor’s series and if 8 is assumed 


* Primed and double-primed coefficients are related to regions 
2and 4, respectively. 


Then to the order of a’, 


Fic. 2. The thick airfoil. 


small compared to a, then it can be shown thatTf 
(neglecting the small quantities in the denominator) : 


p= {2 (2 + DG,)/26, | at + 0(a°) 


where all quantities are evaluated in region 1 ahead of 
the airfoil, From the Prandtl-Meyer differential 
equation for isentropic waves 


dp/d0 = pw*/V M? — 1 
and Bernoulli’s equation 
[(y — 1)/2}w* + a? = a? 
it can be shown that 
dM/dd = —M{[(y — 1)/2])M? + 1}/V 1 


tIn R. & M. No. 1930‘ only the coefficients C; and C; were 
used in deriving an expression for the deflection angle 8. This is 
equivalent to the assumption that all the waves are isentropic, 
and therefore no deflection can be expected. The actual mathe- 
matical error seems to be the omission of the factor '/: from the 
second derivative term of the Taylor expansion of the coefficients. 
Insertion of this factor into the equation derived resulted in a 
zero value of 8 to the order (a*) considered. 


10y Sly + 59 ape 


2 
(—6y? + 407 — 62)M* + (—24y + 60)M? — 24| 


Values of 8/a‘ calculated from the preceding equation for y = 1.40 are given in the following table: 


1.8 


2.0 3.0 4.0 5.0 10 15 20 50 


—332 —10.1 


0.37 1.32 1.37 1.32 1.27 1.31 2.17 3.94 6.75 44.8 146 341 


5,267 


The results indicate a downwash for values of M, less than about 1.28 and an upwash for higher Mach Numbers. 
The magnitudes of the flow deflections between Mach Numbers of 1.2 to 5.0 are extremely small, being always less 


than 0.06° for an a of 0.1 rad. (5.73°). 


and M, = 5, and the result was within 1 per cent of the value given by the approximate theory. 


An exact calculation of the deflection angle was carried out for a = 10° 


It should be 


kept in mind that these results are valid subject to the assumptions that the leading-edge shock is attached and 
that the pressure has not reached a value of absolute zero anywhere in the flow field. 
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The Airfoil With Thickness 


In considering the flow immediately rearward of the thick airfoil, it will be assumed that the configuration js 
such that four oblique shock waves are attached to the nose and tail of the airfoil and that these are the only 
shock waves in the flow field. The flow picture is shown in Fig. 2. Following these assumptions, the flow angle to 


first order can be shown to be 
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— — (& + + [2 ral ~ (a+ + &)* 
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B= 


[(5s + a@)(a@ + + (53 — 6; — + — 


Tue RETURN OF THE .WAKE TO THE FREE-STREAM 
DIRECTION 


Since the forces acting on the airfoil are of finite mag- 
nitude, it is evident from momentum considerations that 
for a two-dimensional flow the downwash far behind 
the airfoil must vanish. Thus the wake returns even- 
The question 
arises as to the mechanism by which the flow just be- 
hind the airfoil is turned back to the free-stream direc- 


It is known that the 
expansion waves in interacting with the shock waves 
neutralize them, until at some finite distance above and 
below the airfoil a complete cancellation has taken 


2C; 
dD 
2C, 
\ 
Nf 
tually to the free-stream direction. 
y tion. 
aoe Consider the flat plate airfoil. 
Fic. 3. System of waves in wake of flat plate airfoil. 
When 


6 = 6 = & = 6; = by 


as for the symmetrical airfoil with equal eo and 
trailing-edge included angles, then 


or 
6 ‘ 6 4 
B= + 1) (: 1) plate 
a Qa 


It appears, then, that the flow deflection fora thick air- 
foil is greater than that of the flat plate at a correspond- 
ing angle of attack and flight Mach Number. For 
example, the angle behind a symmetrical airfoil with 
10° leading- and trailing-edge included angle at a = 3° 
will be about seven times that of the flat plate at that 
angle of attack. 


place. During the process of interaction, however, 
waves are reflected from the points of intersection of 
the rarefactions and the shock (see Fig. 3); also arising 
from these points are lines of entropy discontinuity or 
vortex lines. These reflected waves in turn will be 
transmitted and reflected upon striking entropy dis- 
continuity lines. Somewhere downstream of the air- 
foil, this wave system will have neutralized itself, and 
downstream of that point, the flow will be parallel to the 
free-stream direction. 
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SUMMARY 


The stress distribution in monocoque structures having cut- 
outs has been previously investigated by several authors under 
the assumption of rigid frames. This same problem is here dealt 
with for the circular cylinder reinforced by equal and equally 
spaced stringers and by equal flexible circular frames (rings). 
A solution is developed that considers the deflections of the 
frames according to the theory of curved beams. The results of 
gumerical computations show that the flexibility of the frames 
considerably influences the stress distribution, producing a marked 
stress concentration near the cut. A peak in the shear flows, 
not detected by the previous approximate analyses, appears close 
beside the opening. The test results dbtained by Langhaar and 
Smith’ agree with the more accurate theory within the experi- 
mental scatter limits. 


INTRODUCTION 


- STATE OF STRESS in a monocoque structure 
under any loading that produces shear stresses in 
a given panel of the skin is considerably modified if this 
panel is cut away. The corrective stresses (difference 
between the actual stresses in the structure with cut- 
out and the basic stresses that would occur in the uncut 
skin) may be thought of as perturbations arising from 
the cut and, thence, propagating in the structure. 
These perturbations are here investigated considering 
a factor hitherto neglected which turns out to be of 
paramount importance—namely, the flexibility of the 
frames. 


Review of Previous Investigations 


For the calculation of the stresses in the open por- 
tion of a monocoque structure, some simple graphical 
and numerical methods were worked out by Beskin.! 
The procedure is derived from Wagner’s theory’ of 
torsion bending, further developed by Kappus,’ 
Goodier, von Karman,‘ and others. According to this 
theory the normal stresses in the stringers are propor- 
tional to the area swept by the radius vector from a 
proper pole to the skin midline, the pole being so chosen 
that the normal stresses satisfy the moment equilibrium. 
This procedure is mathematically exact if the shear de- 
formations of the skin and the deflections of the frames 
are disregarded and only the deformations of the string- 
ers in the cutout interval are taken into account. 

Another approximate method was developed by Golo- 
lobov® for monocoque structures with cutout subjected 
to torsion. The shear stresses along the open section 


Received August 28, 1947. 
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Effects of Cutouts in Semimonocoque 
Structures 


P. CICALA* 
Politecnico di Torino, Italy 


are expressed by superposing a constant shear flow on 
the stresses given by the ordinary beam shear theory, 
the two parts being so arranged as to equilibrate the 
applied torque.t It may be shown that for the case of 
the circular cylinder with equal and equally spaced 
stringers, this same solution is obtained by taking 
account only of the shear deflections of the skin in the 
open portion of the structure, all the other elements 
being considered as rigid. 

A more accurate analysis of the effect of cutouts in a 
circular monocoque cylinder was developed by Lang- 
haar and Smith. The deformations of the skin incor- 
porating the stringers are considered for the open bay, 
whereas the end rings of this bay and the structure 
beyond these elements are assumed to be rigid; under 
these assumptions the end sections must remain plane 
and rigid. The basic equation is obtained by mini- 
mizing the strain energy.t By use of this equation 
the authors calculate the stresses in the cutout interval 
and compare the theoretical predictions with test 
results. 

The writer developed® a general procedure for the 
calculation of corrective stresses in semimonocoque 
structures and applied it to the evaluation of the per- 
turbation arising from a cutout in a flat or slightly 
curved wall. Since the same method is applied in 
the present paper, no further details need be given 
here. 


Scope of Present Investigation 


In this paper the effects of cutouts are investigated for 
the circular cylinder reinforced by equal and equally 
spaced stringers and circularrings. In order to separate 
the effect being investigated from other factors, no 
sources of perturbations (e.g., constrained or free sec- 
tions) are supposed to exist in the neighborhood of the 
opening. In other words, the structure extends with 
constant dimensions on both sides of the cut. The skin 
is assumed to be only shear resisting; the stringers, 
axially resisting. The rings are supposed to have no 
stiffness out of their own plane; their flexibility in this 
plane is taken into account according to the theory of 
beams of large curvature. The distance of the centroid 
of ring cross ‘section from the skin cylinder is also con- 


{ The results do not depend on the area of the stringers border- 
ing the cut. This fact is not evident from Gololobov’s analysis, 
which could be considerably simplified. 

¢ This equation of the membrane theory of shells had been pre- 
viously applied to analogous problems.’ 
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Fic. 1. Notations (a), (b), (c) enlarged view showing direction of forces at point P. 


sidered. In this calculation results obtained in pre- 
vious papers*.’ are used. 

The effect of a cutout limited to a single panel or 
extending to many panels is investigated. The influ- 
ence of edge-reinforcements framing the cut is taken 
into account. 


THEORY 


Principal Notations and Hypotheses 


Let z be the axis and R be the radius of the circular 
cylinder covered by the skin, which is supposed to be 
effective in shear with rigidity G and to have constant 
thickness ¢. Let the skin be reinforced by m equal 
and equally spaced stringers of constant area F, axially 
resisting with modulus £, and by equal circular frames 
at constant spacing a. The bay containing the cut 
has length 2a. Let F, be the area of the (constant) 
cross section of the frames, R, the distance of its cen- 
troid G from the z-axis (Fig. 1), and y the coordinate 
measured radially (positive outwards) from the axis 
z, parallel to z through G. Only the circumferential 
stresses in the rings are considered, and plane sections 
(through z) are supposed to remain plane. Therefore, 
these stresses may be expressed by 


My M 
J” RR 
where M, M, denote the moments about the axes z and 


2, Tespectively (positive when giving tension in the 
outer fibers), and 


R+y 


the integration being extended over the area of the 
frame section. The moments are related to the normal 


force N on the same section by 


M, = M — NR, 


Let the stringers be numbered from 7 = 1 to j7 = m- 
Their axial load is denoted by S, positive when ten. 
sion. 

Let the origin of the angular coordinate y¥ (about z) 
be chosen halfway between stringer m and its con- 
secutive 7 = 1. ‘Therefore the coordinate y, of stringer 
jis 

= G — '/2)6 


where 6 = 24/m. 


Let yj’ = (j — l)éandc = Ré. Let the frames be 
numbered from the one at station z = d» (frame 0), so 
that z = ad + as is the coordinate of frame s. At 
station z = —dp lies frame 0’. The panel enclosed 
between stringers 7 — 1, 7 and frames 0, 0’ is termed 
panel 7. The shear flow (tangential stress _ thick- 
ness), constant in each panel, is denoted by g and is 
considered as positive when corresponding to positive 
twist (i.e., rotation of frame s with respect to frame 
s — 1 in the positive y direction). 

Let m be an integer ranging from 2 to (m — 1)/2 
for m odd, to m/2 for m even. By )> is denoted a 

n 


summation extended to the above values of the argu- 
ment » with the convention that, if m is even, the last 
summand corresponding to m = m/2 must be halved.* 


The General Expression of Perturbations for the Semi- 
Infinite Structure 
The structure is here supposed to extend from z = a 
toz = ©. Leta self-equilibrating load system, anti- 
symmetrical with respect to the plane y = 0 through gz, 
be applied to the terminal section z = a. For any 
loading of this type the normal force resisted by 


* This rule avoids a special definition for the values that refer 
to this particular term. The peculiarity derives from the fact 
that the summation of sin? ny, from 7 = 1 to 7 = m has the value 
m/2 for any n except for n = m/2, for which case it takes the 
value m. 
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stringer j at the crossing with frame s may be expressed 


by 

S= (Aint + Aux) sin my; (1) 
the summation being extended as above specified. 
The coefficients A;, Ay; depend upon the applied load- 
ing, Whereas the numbers x1, x11 depend only on the 
properties of the structure. These quantities de-, 
pend also on n; therefore, if confusion can arise, they 
are to be distinguished by the subscript n. 

The numbers x;, x1; may have complex conjugate 
values. In any case their absolute value is <1; this 
means that the perturbation vanishes with increasing 
distance from the loaded section. If x;, x: are not 
real, the coefficients Ay, Ay also become complex 
conjugate. In order to deal with quantities that 
always have real values, it is convenient to introduce 
the quantities 


L = x1 Xm L’ =x + xn : be 
A=A,+An, B = (A; — An)(x1 — xm) 


As shown in previous papers® ® 
=(2+3 1 — L’) 


where* 


4 1 i 
44+-+4+=— -(4+4+-+4+— 
L’ 


(2Ga? t/EFc)(1 — cos né) 


_ Get 2+ cos nd + 
EFa 6(1 — cosns)LR, | JR 


The coefficient f, is a function of and n. 
the approximate expression 


~ n?/(n? — 1) 
may be adopted. More precisely 


tn \= 5 — 6 + (6 — sin 6) 


2 + cos né 
( 1 — cos 6 ) 
cos 6 — cos né 


This equation yields the following values: 


Form = 8: fe = 1.3289; fs = 1.1145; fs = 1.0658. 

For m = 12: fe = 1.3326; fs = 1.12335; fa = 
1.0633; fs = 1.0370; fe = 1.0282. 

For m = 16: fz = 1.3331; fs = 1.1246; fi = 1.0658; 
fs = 1.0401; fe = 1.0262; fr = 1.0184; fe = 1.0157. 

For m = 24: fe = 1.3333; fs = 1.1249; fa = 1.0665; 
fs = 1.0414; fe = 1.0282; fr = 1.0203; fe = 1.0151; 
fo = 1.0115; fio = 1.0090; fi = 1.0074; fiz = 1.0069. 


| 


(Refn — R) | 


In general, 


s the 


* The quantities A, B, L, L’, L”, X, p should be distinguished 
by the subscript , since they all depend on n. 
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In the case of rigid frames, 

L=0, L'= (4) 
1 + (4/3) + VA + (7/12) 

When Eq. (1) applies, the shear flow in the panel en- 


closed between stringers 7 — 1, j and rings s — l, s 
may be expressed by 


q = — xr) + — xn’)] X 


cos ny,’ 
2a sin (né/2) 


(5) 


The state of stress in any ring is completely known 
when the values of moment M and normal force N are 
evaluated for the sections at crossing with stringers. 
For the section of ring s containing the axis of stringer 


Js 
sin ny, 


_ sin ny; 
"1 — cos nd 


M = — 
2a 


(6) 


2a 


n 
cos 6 — cos né 
n 


(6) 


where 


Cy = Ay(2xy? — — + 
Ay — — 


In the intervals between stringers, M varies linearly 
in terms of ¥, whereas the variation of N as a function 
of Wis sinusoidal: viz., for << ¥ < vy 


sin 6 = N(y,) sin — + 
N(¥,-1) sin (¥, — 


The above formulas give the stresses in the whole 
structure, provided the coefficients A corresponding to 
the loading considered are known. 


The Structure Without Opening Subjected to Tangential 

Loading Along the Edge of a Skin Panel 

Let z = 0 be plane of symmetry of the structure 
(Fig. 1) that extends from z = —~ toz = ~. The 
central bay goes from z = —do to z = dp; elsewhere 
the frames have constant spacingt a. The skin is 
supposed to be everywhere continuous. 

Apply along the contour of panel 1 a tangential dis- 
tribution of forces corresponding to a shear flow g* (see 
Fig. 1). The stresses that this loading produces are 
-antisymmetrical with respect to the plane z = 0, as 
well as to the orthogonal plane y = 0. From frame 0 
on they may be represented by Egs. (1), (5), and (6). 
At station z = d» the axial load resisted by stringer 
jis 

+ Since in the cases considered the stresses are antisymmetrical 
with respect to the plane z = 0, the results hold unchanged if a 
frame is located at z = 0. They may be extended to cover ap- 
proximately the case where frames 0, 0’ are stiffer than the 
others by properly increasing the effective spacing a. 
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S, = DA sin ny, (7) 
If g, is the shear flow in panel j, for equilibrium 

— = (8) 


As can be verified by substitution, this equation is 
satisfied by 


= Alcos ny,'/2ao sin (n6/2) ] (9) 
" 

It should be noted that the q value that this equa- 
tion furnishes for panel 1 (yi’ = 0) represents the flow 
on the bordering elements (stringers and rings), whereas 
the flow resisted by the sheet is g, — g* (see Fig. 1a). 

If frame 0, which on its negative side is subjected to 
the flow g, is identical to the other frames, there must 
be the relation* 


ado a XI a XII 
Hence it follows that 


B/A = + 2L(a:/a) (10) 
where ad; = a — dp. 


Furthermore, with somewhat lengthy calculations 
outlined in Appendix I there is derived the equation 


A = 2H,aq* sin (n5/2) (11) 
where 
1 _ my + ala — (11a) 
H, aga(1 — L’ + L)? 6a? 


Eq. (11) that holds true for any value of m (within the 
specified limits) relates the coefficients A to the applied 
qg*. After expressing the coefficients A;, Ay in terms 
of the A, B’s the stresses in any element of the struc- 
ture may be calculated from Eqs. (1), (5), (6), (10), and 
(11). Sothere results 


= g* > H, cos ny,’ (12) 


= 2aq* >> Hy sin (né/2) sin np, (13) 


The shear flow in the panel between stringers j — 1, 7 
and rings 0, 1 is 


* 
= + aol! — a) cos ny,’ (14) 
ain 


The stresses in ring 0 may be obtained from Eq. (6) 
with 
Ca = Ang*(a + do — aL’ — a,L)/a (15) 


If a = 2a the expression of H, simplifies to 


= (4/md)[(1 — L’ + L)*/(1—L)(1 + L’ + L) 


* See reference 9, Eq. 229. This equation states that the shear 
flow given, by Eq. (5) for s = 0 must equal the g value from Eq. 
(9). 
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Fora = a 
= (2/md)[((1 — L’ + — L)] (a7) 


It may be shown tbat similar equations hold if the 
structure outside the central bay is supposed to be 
rigid. If the deflections of rings 0, 0’ are considered, 
we find 


Hy, = 24/m(12 + d + 24p) (18) 


whereas, if these elements are supposed to be rigid, 
HA, = 24/m(12 + d) (19) 
It is shown in Appendix II that the solution expressed 
by Eqs. (12) and (19) might be obtained by a procedure 
similar to that followed in reference 6. The solutions 


indicated in references 1-4 and reference 5 may be de- 
rived as particular cases. 


The Structure with Framed Opening 


The above solution can be used to determine the 
stress distribution in the cylinder near an opening 
whose contour is framed by a stiffening element (edge- 
frame). Suppose that the cut is limited to a single 
panel. The stresses in the structure are first calculated 
under the assumption of continuous skin. In this 
calculation the elementary theory or a more accurate 
one may be used: this makes no difference in the de- 
termination of the perturbation arising from the cut- 
out. Let g be the shear flow so calculated for the 
fictitious panel. This same stress distribution may be 
created in the real structure if, along the framed border 
of the cut, an additional tangential flow, g,, is applied 
of magnitude sufficient to produce the same shear de- 
flection for the edge-frame as for the fictitious sheet 
panel:f therefore, 


(go + = (20) 


where Ky, K, are the shear stiffnesses of the edge-frame 
and of the sheet panel. In a second analysis the 
structure is considered as subjected to the reversed 
imaginary load (Fig. lc). This loading is sustained 
partly by the edge-frame and partly by the surround- 
ing structure. The stresses acting on the bordering 
elements (uncut stringers and rings) are represented by 
the shear flow q: satisfying the equation 


(qr — = (21) 


which equates the deflections of the edge-frame and the 
deflections of the surrounding structure, whose stiffness 
is denoted by K,;. In order to apply the results of the 
previous section to the problem in question, the q 
value from Eq. (21) is to be equated to = value from 
Eq. (12). It follows that 
t See Fig. lb. The stiffening effect of the edge-frame on the 
bordering stringers is here neglected. 
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Hence, g* and, consequently, the perturbation may 
be determined. 

In the particular case where the edge-frame has the 
same stiffness as the replaced panel, g* = O0—that is, 
no perturbation arises. If the cutout is not framed, 
m = —Q- In this case the corrective stresses, when 
added to the basic distribution, give zero stresses for 
the cut panel. 

Particular consideration must be devoted to the 
definition of the stiffnesses K. Consider a rectangular 
wall element (sheet panel or edge-frame) having side 
lengthsaandc. If during the deformation the opposite 
sides remained straight and parallel, the shear deflection 
could be measured by the angle formed by the two 
opposite sides with the normal to the other two. Ina 
general case, where the rectangular element is curved 
on a cylindrical surface and undergoes an arbitrary 
deformation, the shear deflection may be defined by the 
virtual work done during this deformation by a tan- 
gential distribution of forces along the contour cor- 
responding to a shear flow of magnitude 1/ac. By this 
convention, the deformation of a uniformly shear loaded 
sheet panel is equal to the shear strain as normally de- 
fined. The stiffness (ratio of applied flow to corre- 
sponding deflection) for the sheet panel is 


K, = Gt (23) 


For a rectangular plane frame having two sides of 
length a and stiffness EJ’ and two sides of length c and 
flexural stiffness EJ”, 

1/Ky = (a*c/24EJ’) + (ac?/24EJ") 


The calculation of the stiffness of the structure 
surrounding the hole may be readily accomplished by 
use of Eq. (21), considering that g* must replace q, for 
K,=K,. Therefore, 


K,/K, ved H,) — 1 


(24) 


(25) 


Cutouts Extending to Many Panels 


A simple application of the superposition principle 
makes it possible to evaluate the effects of openings 
extending to many panels on the basis of the solution 
presented in the previous paragraphs. The procedure 
may be easily shown by means of an example. 

Consider that an opening of 90° angular extent, with 
edge-frame of negligible stiffness, is cut in a circular 
cylinder having m = 12 stringers, subjected to torque T. 
First the perturbation due to the g* loading is calcu- 
lated. Let p, be the shear flow in panel 7 correspond- 
ing to loading g* = 1 along the contour of panel 1. 
For the particular case under consideration (data as 
specified below) these values are plotted on Fig. 2 at the 
left.* If the cutout extends from stringer m — 1 to 


stringer 2, three load systems gm*, q:*, g2* applied to 


* If deflections are not to be calculated, it suffices to know the 
ratios p;/p; for the evaluation of stresses. The scale of the p 
diagram is then immaterial. 
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| 
he 
| | 
1810° +9'0° 
w=0 
Fic. 2. Corrective shear flows under applied torque T = 


2xR*go, for 30° wide cutout (left figure) and 90° wide cutout 
(right figure). Points from Langhaar-Smith’s tests. 


the contour of panels m, 1, 2 must be simultaneously 
introduced. The resultant stresses are represented 


by superposition 


The unknown quantities gm*, g:*, g.* are then deter- 
mined by putting 


da *a =a 


if 2rR*q is the applied torque. When the unknowns 
are calculated, the shear stresses in the cutout interval 
may be deduced from Eq. (26); similar equations for 
the stresses in any element may be derived from Eqs. 
(12) to (15). The values obtained for the shear 
flows in the open bay of the structure considered are 
represented by the step diagram in Fig. 2 at the 
right. 

If the cut extends to two adjacent panels and these 
are equally stressed, Eqs. (12)—(15) hold true provided 
H,, is replaced by 2H, cos (n6é/2) and the origin of the 
angle y is,taken at the center of the cut. In Eq. (25) 
H,, must be replaced by 2H, cos? (n6/2). 

The procedure may be easily extended to cases where 
the basic values of go are not equal for all the cut panels. 
If an edge-reinforcement of relevant stiffness frames 
the cut, an approximate solution may be derived by 
reducing the g; values according to Eq. (22). 


COMPARISON OF THEORETICAL AND EXPERIMENTAL 
RESULTS 


The above solution was applied to the structure re- 
ferred to in reference 6. The structure tested is a 
circular cylinder of Alclad sheet (¢ = 0.027 in., R = 15 
in.) reinforced by m = 24 stringers of cross-section 
area 0.087 sq.in. and '/,- by 15/s-in. external strap rings, 
spaced a = 14in. Therefore, R, = 15.125 in., F, = 
13/39 sq.in. As the depth of the ring is small compared 
to R, the effect of the curvature may be neglected, 
consequently, J/F, = 1/12 sq.in. Assuming that the 
whole skin area corresponding to each stringer con- 
tributes to the total, F = 0.087 + m 15(0.027/12) = 
0.193 sq.in. It was assumed that G/E = 0.385. 
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Fic. 3. Corrective shear flows for a 30° cutout, according to 
various assumptions, 


With these data, for 2a) = a, from Egs. (3) and (16) 


H2= 0.0002s, H; = 0.0018, Hy, = 0.0053 
H,; = 0.0108, Hy, = 0.0161, Hy = 0.0192 (27) 
Hs = 0.0203, Hy = 0.0202, Hy = 0.0198 
Ay, = 0.0195, Hy = 0.0194 

For do = a, from Eq. (17), 
Hy, = 0.0004,, Hs = 0.0025, Hy, = 0.0059 
Hs = Hes = 0.0107, Hz = 0.0109 (28) 
Hs = 0.0105, Hy, = 0.0098, Fy = 0.094 
Hy = 0.0091, Hy = 0.0090 


For 2a) = a, if the frames are assumed to be rigid, 
from Eggs. (4) and (16), 


FH, = 0.05565, Hs = 0.0463, IH, = 0.0391 
Hs = 0.0335, Hs = 0.0293, Hy = 0.0261 (29) 
Hs = 0.0237, H, = 0.0220, Hy = 0.0209 
= 0.0202, ITs = 0.0200 
With Langhaar’s assumptions, from Eq. (19), 
Hz = 0.0786, HH; = 0.0736, H, = 0.0680 
Hs; = 0.0625, IT, = 0.0575, IT; = 0.0532 (30) 


Hs, = 0.0498, = 0.0472, 
Ay = 0.0443, FT = 0.0439 


Golobov’s assumptions give, for any n, H = 0.0833. 

Comparison of values in Eq. (27) with Eq. (28) 
shows that the effect of the flexibility of the frames 
considerably reduces the amplitudes of the lowest har- 
monic components, whereas the higher components re- 
main practically unchanged. On the contrary, the 
effect of the warping of the terminal sections is more 
marked for the higher harmonics: this fact is shown by 
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comparison of values in Eq. (29) with values in Eq. (30), 
which were calculated under Langhaar’s assumption 
that warping was prevented. Comparison of values in 
Eq. (29) with Eq. (28) indicates that, in the case of a 
longer cutout, the frame flexibility has a less marked 
influence on the lower harmonics and a more marked 
influence on the higher harmonics than in the case of a 
shorter cut. 


The shear distribution in the cutout interval, in the 
case of a 14 in. long, 30° wide opening was calculated 
according to the various assumptions and plotted in 
Fig. 3. The scale of the diagrams is defined by their 
common ordinate at the cut. For this case the ratio 
K,/K,, when calculated according to Langhaar’s as- 
sumptions, comes out 1.48. If the warping of the end 
sections is considered (infinitely long structure, rigid 
frames) there results K,/K, = 0.55. If the flexibility 
of the frames is considered, one obtains K,/K, = 
0.081; the flexibility of the structure is considerably 
increased in this case. 


Fig. 4 affords a comparison of the theoretical results 
with test data obtained by Langhaar and Smith in tor- 
sion tests on the specified structure. The test results 
obtained in a bay of the cylinder with uncut skin are 
compared in Fig. 4a, with the theoretical flow value 
go = T/2xR? (T = torque). In Fig. 4b, for a 30° 
by 14-in. cutout, the corrective stresses computed for 
the open bay according to values in Eq. (27) (full line) 
and Eq. (29) (dotted line) are compared with test data. 
These confirm the peak value near the cut and the de- 
creasing value at distance. The points plotted in Fig. 
4c refer to the bay adjacent to the same opening as 
above; the same scale is used as for Fig. 4b. The full 
line is derived from Eq. (14). 


Figs. 4d and 4e represent the corrective shear flows 
in the open bay and in the adjacent one, respectively, 
for a 60° by 14-in. cutout. Full lines are calculated 
according to values in Eq. (27); the dotted line is ob- 
tained with Langhaar’s assumptions. The same scale 
applies to both diagrams, as defined by the value at the 
cut. The test points show qualitative agreement with 
theory. 

The calculations were also performed for the same 
structure, considering the cross-section area to be con- 
centrated in twelve stringers (m = 12, F = 0.772 sq.in.). 
The corrective shear flows so calculated for the bay 
containing a 60° cutout are plotted on Fig. 5. Line 1 
refers to a 14-in. long opening; line 2, to a 28 in. open- 
ing. Line 3 was computed assuming the ring material 
to have a modulus E three times greater (ppew = 1/spota): 
The three lines do not show a marked difference from 
each other, whereas they differ considerably from the 
dotted line obtained according to Langhaar’s assump- 
tions. For the same structure, in the case of a 30° 
by 14-in. cut, the corrective shear flows are plotted in 
Fig. 2 at the left. If the cut extends to 90°, in the case 
of constant basic shear flow gq (for instance, under 
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applied torque 7), the corrective stresses are repre- 
sented by Fig. 2 at the right. The points plotted there 
represent test values from reference 6. Some tests 
were conducted on a somewhat modified structure 
(bordering stringers and rings stiffened, 28-in. instead 
of 14-in. cutout length). Since the theoretical line 
does not account for these changes, only an approxi- 
mate comparison is offered. Nevertheless, it may be 
stated that the fundamental results of the theory 
(namely, peak value near the cut,* vanishing value at 
adistance) are confirmed by the test. 

The same structure was also tested under bending 
and shearing loads. In this case only the stresses in the 


* Some tests indicate a moderate stress in the proximity of the 
cut, contrary to theoretical predictions. This result might be 
ascribed to the edge stringers being not fully effective hathinesed 
as a consequence of local buckling). 
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Fic. 4. Comparison of test results by Langhaar-Smith, with 
theoretical values for torque loading T = 27R%q. (a) Shear 
flows in a bay of the structure without opening. (b) and (c) 
Corrective shear flows in the cutout interval and in the adjacent 
bay, respectively, for 30° wide opening. (d) and (e) Corrective 
shear flows in said bays, for 60° wide opening. 
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Fic. 5. Corrective shear flows for the structure having m = 
12 stringers, 60° cutout. Line 1 refers to 14-in. long opening, 
dural rings; line 2 to 28-in. long opening, dural rings; line 3 
to 14-in. long opening, steel rings. Dotted line refers to Lang- 
haar’s assumptions. 


stringers were measured. Since the corrective stresses 
here are added to high basic stresses, the scatter of the 
test points renders the comparison with theoretical 
predictions less conclusive. 


CONCLUSIONS 


The procedure developed in this paper allows the 
calculation of the corrective stresses representing the 
effect of cutouts in semimonocoque structures of circu- 
lar section. The starting point of the procedure is the 
evaluation of the “unit” perturbation produced by a 
tangential self-equilibrating loading applied along the 
contour of a single panel. This perturbation is repre- 
sented as a summation of harmonic components; the 
stresses in the elements where the highest values occur 
may be computed by use of Eqs. (12)—(15), with the 
aid of Eq. (lla), which in the most frequent cases 
takes the simpler form of Eq. (16) or (17). This per- 
turbation must be added to the basic system of stresses 
which is determined beforehand assuming the skin to 
be continuous. The magnitude of the corrective 
stresses, expressed by the factor g* is derived from Eq. 
(22) when the cut extends to a single panel. When 
many adjacent panels are cut, the method of super- 
position permits the computation of the corrective 
stresses from unit perturbations, as Eq. (26) indicates 
for the case of a three-panel cut. Therefore, the cal- 
culation of the perturbation need not be repeated for 
each loading condition considered. 

The procedire refers to states of stress antisymme- 
trical with respect to the mid-section of the cut. The 
stiffening effect of edge-reinforcements framing the 
cutout is considered, with some approximate assump- 
tions regarding the effect of the reinforcement on the 
behavior of the stringers. Further investigations are 
to be devoted to symmetrical loadings (see asterisk 
footnote on page 178). 

From the results of numerical applications by the 
procedure given, it may be concluded that the as- 
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sumption of rigid frames accepted by previous investi- 
gators may lead to considerable error. The peak values 
in the stress field (especially in the shear flow) that the 
more accurate theory indicates for the elements in the 
neighborhood of the cut reach values many times 


’ greater than the maxima indicated by rigid frame 


theories. 

The problem needs more experimental evidence. 
Langhaar and Smith’s tests show good agreement on the 
fundamental results of the theory. But the experi- 
mental investigation should be pursued in order to as- 
certain whether, when the scatter limits became 
smaller, the agreement with theory would be im- 
proved. 


APPENDIX I 


In order to derive Eq. (11), the procedure here 
sketched may be followed: ' 

Cut the central bay with its terminal frames 0, 0’ 

. free from the structure. 

Apply at frame 0 the axial forces S; = sin my, to the 
stringers. 

Calling o,, 7, the normal and tangential stresses due 
to this imaginary loading, calculate the virtual work 
done by these forces during the deformation produced 
by the g* loading, with the aid of 


= w, sin ny, 


j= 


where the first integral is taken over the volume of 
stringers, stressed according to Eq. (7), and over the 
volume of rings 0, 0’ where the o’s may be computed 
according to Eqs. (6). The second integral is taken 
over the sheet volume, where the r’s may be derived 
from Eq. (9); w,; is the axial displacement of the loaded 
section of stringer 7 expressed by 


—2EFw, =a (4 + Ay =) sin ny; 


according to Eq. [106].° 

Taking account of the orthogonal properties of sine 
functions and the above equations, Eq. (11) may be 
derived. 


APPENDIX II 


If only the deformations of stringers and sheet form- 
ing the open bay are taken into account, the terminal 
sections are so constrained that they remain plane and 
rigid (Langhaar’s ‘assumptions). Assume first the 
skin continuous and apply the g* forces as in the pre- 
ceding analysis (see Fig. 1). The corresponding stresses 


* Under the assumption of rigid frames, the effect of openings 
in plane or curved reinforced walls (stress diffusion problem) has 
been investigated for symmetrical loadings (reference 9, para- 
graphs 19 and 20). 


satisfy the following equationT 
(Ga*t/6EFc)(S; — S;-1) + gy = (31) 


where a = 2a) is the length of the bay and, for j = 


v’ + v” cos 

v’ + v" + g* 

with v’, v” constants (proportional to the relative rota- 
tion and displacement at center of the terminal sec. 
tions). Considering that the forces S and the flows g 


form self-equilibrating systems, from Eq. (31) it follows 
that 


0; 


™m™ ™m 
v= v, cos = 0 (32) 
Consequently, it may be written 


vy = (29*/m) 2 cos ny,’ (33) 


Substituting for S, g, v from Eqs. (7), (9) and (33) 
into Eq. (31) and equating the coefficients of cos 
ny,’ on both sides, we obtain 


A [1+ (4/12)] = (2/m) ag* sin (n5/2) 
in accordance with Eqs. (11) and (19). 

The solution of this same problem may be given an- 
other form. Substitution of Eq. (8) into Eq. (31) 
yields 

(Ga*t/12EFc) (2g; — + = (34) 

The solution of this equation may be written 

= C(x, + x™4%) + C’ + cosy,’ (34a) 
where x satisfies the equation 


x + (1/x) —2 = 12EFc/Ga't 


and the constants C; C,’ C” are determined from the 
known value of g; and the equilibrium conditions 


= cos y,’ = 0 (35) 


For Ga*t/EFc = 0 (Gololobov’s assumptions), the 
solution takes the form 


a= C+C'+C", = C’ + C" cos yy’ 


(j = 2, 3,...,m) 


For EFc/Ga*t = 0 (Wagner’s assumptions), the solu- 
tion may be put in the form 


In any case Eqs. (35) apply. 
If the stringers m and 1 have area F’ # F, Eq. (34) 
for j = 2 (andj = m) is modified to become 


+ This equation may be derived from Eq. (88)* by considering 
that the same constraint conditions as in the real structure would 
be attained if the end sections of the bay were planes of sym- 
metry. 
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CUTOUTS IN -SEMIMONOCOQUE STRUCTURES 


(36) 


( F’ F 

The solution of Eq. (34a) still holds; also Eqs. (35) 
apply. Eq. (36) substitutes the third condition de- 
fining the constants. 
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Letters to 


Dear Sir: 

In my paper ‘‘Numerical Procedures for Calculation of Stress 
Distribution in Sheet-Stringer-Frame Combinations’? (JOURNAL 
OF THE AERONAUTICAL SCIENCE, Vol. 14, No. 8, p. 481, August, 
1947) no startling or new discoveries were intended. Use was 
made of the most fundamental concepts in basic elasticity. The 
standard method of finite differénces was employed to reduce a 
linear differential equation to a determined set of linear algebraic 
equations. The method chosen for solution of these equations 
was essentially that of matrix iteration. For physical clarity, 
the operational form of these calculations was made identical 
with that of the Shortley, Weller, and Fried article of 1940. 

What then was the purpose of an article of this type? 

Although there are many numerical methods available for 
shear lag analysis, it appeared that several of the major ones were 
identical in their basic assumptions and led to recursion or re- 
laxation formulas that were reducible to a common form: e.g., 
Hoff (PIBAL), Kuhn (N.A.C.A.), Kempner (N.A.C.A.), Wignot 
(Lockheed), etc. It was felt that perhaps this derivation, being 
more fundamental and more clearly motivated (mathematically), 
brought out the basic similarity among the methods. Further, 
it allowed direct extension on a common basis to the problems of 
bending and torsion of a box beam, cutouts, and swept box beams 
(as yet unpublished). 

In the simplest case of shear lag, all of these methods neces- 
sarily involve somewhat similar dimensionless parameters. To 
my knowledge, J. Wignot’s report first incorporated the use of 
the term ‘‘characteristic length’ and first referred to this di- 
mensionless parameter as §?. (In Kuhn’s notation this would be 
K/1?.) 

I regret the oversight on my part in not including references for 
these terms. In view of the fact that it is intended to use this 
parameter and this ‘‘concept’”’ for comparing structures on a non- 
dimensional basis, I should like to emphasize that any credit for 
doing this belongs solely to the structural methods group at 
Lockheed. My interests are confined to analytical procedures, 
and the sole purpose of the article was to show the application of 
standard procedures to a standard type of boundary value prob- 
lem in mathematical physics. I believe that this application 
has never before appeared in the literature. 


MARVIN STERN 
Structures Engineer 
Republic Aviation Corporation 


the Editor 


Dear Sir: 

This writer wishes to point out an interesting analogy that 
exists in the calculation of the intersection of characteristics in 
two-dimensional steady isentropic flow and one-dimensional 
unsteady isentropic flow. The solution in the steady flow has re- 
cently been described by Cronvich' but will be repeated here to 
make the analogy clear. Consider two characteristics, A and B, 
in a steady flow which intersect, and are refracted as, waves A’ 
and B’ (see Fig. 1). The rule for the calculation of the condi- 
tions at region 4 is that the change of the angular quantity v across 
an incident characteristic is equal to the change of v across the 
refracted characteristic, or 


and 

The angle v is the angle that a two-dimensional steady isen- 
tropic flow must be deflected in an expansion around a bend to 
reach a Mach Number, M, when starting initially from a Mach 
Number of 1.0. The theory tells us that » is only a function of 
Mach Number, this function being 


y¥+1 
= ta M? — 1) — 
arc tan 1M? — 1 


Fic. 1. 


| 

34a) 
7 
35) 
the 
m) 
4) B 3 A 
ng 
m- 
° 


JOURNAL OF THE AERONAUTICAL SCIENCES+~MARCH, 


Fic. 2. 


Also, since the change in v across a characteristic is equal to the 
angle of the flow deflection in passing through the wave, it can be 
seen that with this method the solution of such interactions is a 
simple matter. It should be pointed out that the above method 
of solution was worked out independently by A. R. Kantrowitz 
in an unpublished paper. 

Consider now the intersection of two characteristics in a one- 
dimensional unsteady isentropic flow, as depicted on the time- 
distance, ¢-x plane, of Fig. 2. As before, the incident character- 
istics are A and B, and the refracted characteristics are A’ and 
B’. The Riemannian? theory of wave propagation tells us that, 
for this type of flow, in going across the upstream moving 
wave BB’, the parameter P remains constant, and in going 
across the downstream moving wave AA’, Q remains constant, 
where 


P=u+t [2/(y —1)]a 
Q=u — [2/(y — 1)Ja 


Thus, to solve for the values of velocity u and speed of sougd a 
in the region 4, we write 


1948 


us + 


and 


us — = — a 

3 1” 4 

where conditions in regions 2 and 3 are already known. These 
equations can easily be solved to obtain the desired result, 
However, if the above equations are.solved using also the fol. 


lowing 


and 


it can easily be shown that a rule analogous to that used in the 
steady intersection problem may be formulated as follows: 
The change of velocity and speed of sound across the incident leg of g 
characteristic 1s equal to the change of velocity and speed of sound 
respectively across the refracted leg. This rule may be stated 
mathematically as follows: 


U2 = 

a a2 a3—- 
and 

UW U3 = 

a = 
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Note on the Simple Ram-Air Intake 


Preceded by anu Shock in Supersonic 
Flight 


‘ JOHN D. STANITZ* 


National Advisory Committee for Aeronautics 


SUMMARY 


This note considers the performance of a simple ram-air 
jntake at supersonic flight speeds. The intake is preceded by a 
normal shock after which the flow is diffused at subsonic veloc- 
ities. Exact expressions are developed for the ram pressure 
ratio and the ram efficiency as functions of the flight Mach 
Number and the diffuser efficiency. In addition, a simple, 
approximate equation is developed for the ram pressure ratio 
at flight Mach Numbers greater than 2.0. The conclusions of 
this analysis are: (1) For flight Mach Numbers greater than 
2.0, the ram pressure ratio varies approximately with the square 
of the Mach Number; (2) for a diffuser efficiency of 1.0, the 
minimum ram efficiency is 0.76; (3) for very high flight Mach 
Numbers, the ram efficiency becomes asymptotic to 1.0; and 
(4) for diffuser efficiencies less than 1.0, the ram efficiency is 
improved by normal shock compression at low flight Mach 
Numbers. 


SYMBOLS 


inlet area of ram-air intake 

Mach Number relative to ram-air intake 
static pressure 

flight velocity 

flow rate through ram-air intake 

ratio of specific heats 

ram efficiency 

diffuser efficiency 

density of atmospheric air 


Subscripts 


0, 1, 2 refer to stations shown in Fig. 2 
t refers to total pressure 


INTRODUCTION 


Dyers THE FLOW OF AIR through a simple ram- 
air intake at supersonic flight speeds. Let the 
inlet area, A, be simply connected and normal to the 
path of flight. If the required flow rate is 


W = pAV (1) 


where p is the atmospheric density and V is the flight 
velocity, the air enters the intake in a manner similar 
to that shown in Fig. 1. An oblique shock originates 
externally at the leading edge. The internal flow, 
beyond the initial disturbance lines shown in Fig. 1, 
depends upon the configuration of the intake, which 


need not be considered in this analysis. 


Received September 24, 1947. 
* Mechanical Engineer, Turbine and Compressor Division, 
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If the required flow rate is less than the amount 
specified by Eq. (1), then from continuity considera- 
tions it is necessary for the excess air to be deflected 
around the leading edge of the intake. However, 
from the rule of forbidden signals,’ supersanjc flow 
cannot be deflected in this manner, and it ‘becomes 
necessary for the oblique shock to detach ite fom 
the leading edge of the intake and to form an essentially 
normal shock ahead of the intake so that the flow can 
be deflected at subsonic velocities as shown insFig, 2. 


Fic. 1. Flow: for intake with air 


flow rate equal to p. AV. 


Fic. 2. Flow configuration for simple ram-air intake with ‘aie 
flow rate less than pA V. 
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Fic. 3. Ram pressure ratio, pi2/po, as a function of the flight 
Mach number, Mo, for several values of the diffuser efficiency, 
np. Ratio of specific heats 1.4. 


This normal shock ahead of the simple ram-air intake 
occurs no matter what the internal configuration of the 
intake, provided the required flow rate is less than that 
specified by Eq. (1). In this note the performance of 
a simple ram-air intake preceded by a normal shock 
at supersonic speeds has been investigated. 


ANALYSIS 


The ram efficiency, nz, is defined as the ratio of the 
kinetic energy of the air leaving the ram-air intake, 
if expanded isentropically to the free-stream pressure, 
to the kinetic energy of the air approaching the intake. 
Expressed in terms of the static and total pressures 
this definition becomes 


me = [1 — — (p0/p (2a) 


where is the static pressure, 7 is the ratio of specific 
heats, the subscript / refers to total pressure, and the 
numerical subscripts refer to the stations shown in Fig. 
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The denominator of Eq. (2a) can be expressed jn 
terms of the flight Mach Number so that 


/ (po/Pa) 
Mo*/{[2/(y — 1)] + 


The ram pressure ratio, pa/po, is expressed by the 


identity 
Po Pi /\Pu/\P0/\bo 


and is therefore equal to the product of four pressure 
ratios, each of which can be determined. 

The diffuser pressure ratio, p2/pi, is given by the 
equation 


= {11+ wl(y — 


where np is the diffuser efficiency of the ram-air intake 
as defined by Eq. (4). 
The ratio of the total to static pressure at station 1 is 
given by 
= (1 + — Ga) 
and at station 0 
= {1 + [(y — Gb) 


The ratio of total pressures across the normal shock, 
Pu/Pw, depends only upon the flight Mach Number, 
Mb, and is given by Eq. (6):? 


2 ty) 
Po \y+1 
— 1) 6) 
(y + 1) Mo? 
The normal shock ahead of ‘the ram-air intake 
reduces the air velocity from supersonic ahead of the 


NR (2b) 


shock (station 0) to subsonic after the shock (station . 


1). This change in Mach Number is given by Eq. (7): 


— [(y — 1)/2] 


The ram pressure ratio, p»/po, is now obtained by 


combining Eqs. (3), (4), (5a), (5b), (6), and (7), 
reducing to 


Po ¥+1 


2 


{[27/(y — — 


(7) 


- (8) 


This equation determines the ram pressure ratio in 
terms of the flight Mach Number, Mo, and the diffuser 
efficiency, np. 

The ram efficiency, nz, is obtained directly from 
Eqs. (2b) and (8): 
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The ram pressure ratio, p2/po, given by Eg. (8), is 
shown by the solid curves in Fig. 3 for y equal to 1.4 
and for diffuser efficiencies equal to 1.0, 0.6, and 0. 
As the flight Mach Number increases, the slope of 
these curves approaches 2.0, so that the ram pressure 
ratio increases approximately with the square of the 
flight Mach Number. This fact is used to obtain a 
simple expression for the ram pressure ratio at flight 
Mach Numbers above 2.0. 

First consider the case for a diffuser efficiency 
equal to zero. Eq, (8), with y equal to 1.4, then reduces 


to 


RAM PRESSURE RATIO 


bo = — 1)/6 

and since, for Mo greater than 2.0, 
7M,? >> | 

the ram pressure ratio becomes 

Pa/po ~ (7/6)Mo? 


For a diffuser efficiency of 1.0 the ram pressure ratio 
is more nearly approximated by the empirical equa- 
tion 


(10) 


P2/ po 1.33M_? (8/6) M,? 


which by comparison with Eq. (10) suggests the 
following approximate equation: 


Po nv) /6\M.? (11) 


This equation is given by the dashed lines in Fig. 3 
for np equal to 1.0, 0.6, and 0. For flight Mach 
Numbers greater than 2.0, Eq. (11) is seen to agree 
closely with the exact values of the ram pressure ratio 
given by Eq. (8). 


Ram EFFICIENCY 


The ram efficiency, ng, given by Eq. (9), is plotted 
in Fig. 4 for y equal to 1.4 and for diffuser efficiencies, 
np, of 1.0, 08, 0.6, and 0. For a diffuser efficiency of 
1.0, the ram efficiency decreases with increasing Mach 
Number until a minimum value of 0.76 is reached at 
a Mach Number of about 5.0, after which the efficiency 
increases. This increase in efficiency at high flight 


Mach Numbers occurs because the normal shock losses 
do not increase so rapidly, with increasing Mach 
Number, as does the energy input. 


an in- 
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Fic. 4. Ram efficiency, 72, as a function of the flight Mach 
number, Mo, for several values of the diffuser saaveat ”pD. 
Ratio of specific heats 1.4. 


spection of Eq. (2a), or Eq. (2b), both the numerator 
and the denominator approach 1.0 as the flight Mach 
Number approaches infinity, so that the efficiency 
curves in Fig. 4 are asymptotic to a value of 1.0 at 
very high flight Mach Numbers. 

For diffuser efficiencies less than 1.0, the ram effi- 
ciency at first increases with flight Mach Number 
until a peak is reached, after which the efficiency 
decreases to a minimum and then rises again, becoming 
asymptotic to 1.0 for very high flight Mach Numbers. 
The initial peak efficiencies occur because, for low 
Mach Numbers, the normal shock efficiency is higher 
than the diffuser efficiency. The normal shock 
efficiency is, in fact, given by the ram efficiency, 
nr, for np equal to 1.0. ; 


CONCLUSIONS 


(1) For flight Mach Numbers greater than 2.0, 
the ram pressure ratio varies approximately with’ the 
square of the flight Mach Number. 

(2) For a diffuser efficiency of 1.0, the minimum 
ram efficiency is 0.76, which occurs at a flight Mach 
Number of about 5.0. 

(3) For diffuser efficiencies less than 1.0, the ram 
efficiency is improved by normal shock compression 
at low flight Mach Numbers. 

(4) For any diffuser efficiency, the ram efficiency 
becomes asymptotic to 1.0 at very high gate Mach 
Numbers. 
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' technical book they wish to borrow. Through an exchange 
‘agreement with the Engifecring Societies Library, any 
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Technical Information Service 
This service has experienced personnel under the super- } 


vision of trained aeronautical engineers to compile any in- 
formation desired. The services range from listing special- 
ized reference books to the preparation of exhaustive 
bibliographies, digesting of reports, and general surveys of } 
any acronautical subject. Some of the available services 
are: 


Bibliographies on any aeronautical subject. 

Reports on any aecronautica] subject. 

Digests of acronautical books, f Papers, periodicals, and refer- 
ences. 

Translations. 

Engineering investigations of special aeronautical subjects. 

Biographies of individuals engaged in aeronautics. 

Photostats of any aeronautical or general engineering 
material. 

Microfilms made on special order. 

Photographs made from thé Institute’s photographic collec- 
tion. 

Drawings and tracings made. 


In, addition to the services mentioned any commission 
which comes within the scope of the Service will be ac- 
cepted. Special arrangements may be made for work re- | 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any material in the collections of 
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